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Uber einige Methoden der Relaxationsrechnung 


Von EDUARD STIEFEL, ETH., Ziirich 


Die von SouTHWELL!) und seinen Mitarbeitern als Verfeinerung des Seidel- 
schen Verfahrens und der Methoden von H. Cross?) entwickelte Relaxations- 
rechnung’) zur iterativen Auflésung von linearen Gleichungen hat sich haupt- 
sichlich fiir die Lésung von Randwertaufgaben bei den partiellen Differenzen- 
gleichungen der mathematischen Physik bewahrt. In der vorliegenden Arbeit 
werden zunachst im Abschnitt 1 die Grundlagen dieser Rechentechnik erneut 
auseinandergesetzt und zugleich die mathematischen Hilfsmittel bereitgestellt, 
die dann in Abschnitt 2 erlauben, das allgemeine Prinzip der Relaxationsrech- 
nung mathematisch zu begriinden. Sie fiihren in Abschnitt 3 zur allgemeinen 
Methode der simultanen Relaxation. 

Spater wird dann dem Rechner die in der Relaxationsrechnung mit Recht 
so geschitzte Freiheit in der Wahl der einzelnen Schritte genommen, und es 
werden —im Hinblick auf den Einsatz von Rechenautomaten — einige zwangs- 
laufige Methoden geschildert. Nach den bekannten Verfahren des starksten 
Abstiegs und der Iteration in Gesamtschritten wird in Abschnitt 5 ein 
«n-Schritt-Verfahren» angegeben, welches sowohl eine sukzessive Approxima- 
tion an die Lésung als auch deren exakte Bestimmung in endlich vielen Schrit- 
ten liefert. Besondere Aufmerksamkeit muBte der Tatsache geschenkt werden, 
daB alle Relaxationsprozesse dazu neigen, am SchluB die Konvergenz durch 
Bildung von «Kafigen» zu verlangsamen*). 

Herrn Dr. F. Krantz sei gedankt fiir seine Entwicklung der numerischen 
Rechentechnik fiir die Methoden der Abschnitte 1 und 2. Seine durchgerech- 
neten Beispiele haben mich oft unterstiitzt. Herrn Dr. H. RuTIsHAUSER und 
Herrn U. HocustrassEr danke ich fiir die Durchfiihrung einiger gréBerer Bei- 
spiele auf der programmegesteuerten Zuse-Rechenmaschine in Ziirich. Herr 


1) R. V. SOUTHWELL, Relaxation Methods in Engineering Science (Clarendon, Oxford 1946) ; 
Relaxation Methods in Theoretical Physics (Clarendon, Oxford 1946). 

2) H. Cross, Analysis of Continuous Frames by Distributing F ixed-End Moments, Proc. Amer. 

Soc. chem. Eng. 1930); Numerical Methods of Analysis in Engineering (Symposium Illinois _In- 
stitute of Technology, Macmillan, New York 1949). 

3) Man vergleiche auch: L. Coirvatz, Numertsche Behandlung von Differentialgleichungen 
(Springer, Berlin 1951), S. 106 und 295. Ferner: G. Tempie, The general Theory of Relaxation 
Methods Applied to Linear Systems, Proc. Roy. Soc. London [A], 169, 476-500 (1989). Dort sind 
das allgemeine Verfahren der Relaxationsrechnung (Abschnitt 2) und das Verfahren des starksten 
Abstieges (Abschnitt 4) geschildert und Konvergenzbeweise gegeben. 

4) Allgemeine Konvergenz- und Fehlertheorien findet man bei: J. VON NeumANN und H. H. 
Go.pstine, Numerical Inverting of Matrices of High Order, Bull. Amer. Math. Soc. 43, Nr. 11 (1947). 
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A EpuUARD STIEFEL ZAMP 


o% - a be 
Dr. E. BASCHELET, Privatdozent in Basel, war so freundlich, mich auf einige 
Literaturstellen aufmerksam zu machen. 


1. Grundlagen 


Zur Bestimmung der 7 unbekannten GréBen 1, us, ..., ¢, Sei das symme- 
trische lineare Gleichungssystem vorgelegt 


D> Vix ty + 1; =0 (i= 1, 2p.e5 %) eee 


k=1 


din = Tei » (2) 


mit den Koeffizienten d,, und den absoluten Gliedern /;. Setzt man in die linken 
Seiten der Gleichungen irgendwelche Naherungswerte fiir die Unbekannten 
ein — wir nennen sie v,, Vs, -.., Un ~, So werden die rechten Seiten Werte 74, 7, 
..., %, annehmen, die im allgemeinen von Null verschieden sind und Restduen 
heiBen: 


n 


D> Fxtrth=n.- (3) 


k=1 


Der Grundgedanke der Relaxationsrechnung besteht nun darin, durch syste- 
matische Korrekturen Av,, Avg, ..., Av, der Versuchswerte die Residuenschritt- 
weise immer kleiner zu machen. Diese Korrekturen ergeben die neuen Residuen 


r, = Ts a SY: din y lv; . (4) 
k=l 


Zur Vereinfachung der Schreibweise fassen wir die Unbekannten und die Nihe- 
rungswerte zu Vektoren “= (1, Mo, ..., %,) und v = (v4, Ug, ..., ¥») ZUSaMMen 
und fiihren auch den Konstantenvektor / = (I,, J, ..., /,) und den Residuen- 
vektor 7 = (7,;) ein. Einem beliebig gegebenen Vektor x = (x, %», ..., %,) kann 
durch die Formeln 


n 


he ay din, Xx (5) 


k=l 


ein neuer Vektor y = (4, Vo, ..., Vn) Zugeordnet werden. Diese Operation 
(lineare Transformation) bezeichnen wir durch das Symbol D, schreiben also 


Bele (6) 


und nennen D den Operator des gegebenen Problems (1). Die Beziehungen (1), 
(3), (4) kénnen damit geschrieben werden 
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Dut 0;, (1’) 
Do pl=7, (3’) 
y=r-+ DAv. (4’) 
Versteht man unter dem skalaren Produkt zweier Vektoren a = (Gin dona 


und} = (b,, b,, ..., b,,) den Ausdruck 


Ww 


(a, b) cS DEG, by , 


R=1 


so ergibt sich fiir das skalare Produkt aus einem Vektor x und seinem Trans- 
formierten y = Dx unter Berticksichtigung von (5) 


(x, Dx) = = 3% v= De DLE ee | 
(7) 
(9G, LD) =) eo he | 


(4, k) 


also die mit den Koeffizienten d,, gebildete quadratische Form. Wir nennen sie 
Q(x) und setzen im folgenden voraus, daB8 sie positiv definit sei, das heiBt 


Q(x) = (%, Dx) > 0 (8) 


fiir jeden Vektor x = (4%, %, .--, %), der vom Nullvektor verschieden ist. 
Ist z ein weiterer Vektor, so findet man noch 


(gD ae iia, (D2 es De din tn = Bo din Hee 


(2) (2, k) 


und andererseits 


(Dx, z) = (Dates Pe ae 
(1) (4, R) 


() (R) 
Wegen (2) folgt daraus die wichtige Beziehung 
(Dx, 2) = (x, Dz) , (9) 


womit auch die Symmetriebedingung (2) symbolisch geschrieben ist. Wir 
bezeichnen (9) auch als Greensche Formel, da sie fiir einen Differenzenoperator 
D das Analogon zur bekannten Greenschen Formel der Differentialrechnung 
darstellt. 

Wir wollen gleich einen solchen Fall zur Illustration behandeln. Es handle 
sich um die Divichletsche Randwertaufgabe der Differenzenrechnung fiir ein 
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gegebenes Grundgebiet, das in Figur 1 quadratisch angenommen wurde. Dieses 
Quadrat wurde schachbrettartig unterteilt, so daB ein quadratisches Gitter 
entsteht. Ist in jedem Gitterpunkt eine Zahl 
angeschrieben, so sprechen wir von einer Gztter- 
junktion. Es sei nun eine solche Gitterfunktion 
u gesucht, so daB 
1. ihr Wert in jedem inneren Gitterpunkt 
das arithmetische Mittel der vier be- 
nachbarten Werte ist und 
2. die Funktion in den Gitterpunkten des 
Randes gegebene Werte annimmt. 
Numeriert man in Figur 1 die inneren 
Gitterpunkte in der Reihenfolge, wie man die 
Zeilen eines Buches liest, so hat man total 25 
Fig 1 Unbekannte 1, Ws, ..., Ye; und erhalt zu ihrer 
Bestimmung 25 Gleichungen. Zum Gitter- 
punkt Nr. 9 gehdrt zum Beispiel die Gleichung 


Ug + Ug t Uyo + Uy, 
2 ‘14 


Ug ie 
oder 
4 Ug — Ug — Ug — Uy — Uy = 0, (10) 


wahrend der Gitterpunkt Nr. 10 zur Gleichung AnlaB gibt 
4 Uy9 — Uy — Uy, — Uz — frp =O. (11) 


Dabei bedeutet fp den gegebenen Randwert im Randpunkte R. Diese 25 Glei- 
chungen haben die Symmetrieeigenschaft (2), so daB alle bisherigen Uberle- 
gungen anwendbar sind. Was wir oben einen Vektor nannten, ist nun hier eine 
Gitterfunktion, genauer eine Werteverteilung auf den inneren Gitterpunkten. 
Die Ausiibung des Operators D auf eine Gitterfunktion x ergibt nun als Wert 
von y = Dx 
im Punkt Nr. 9: 

Vo = 4 Xy — X4 — Xy — Hy — XH, (12) 
im Punkt Nr. 10: 


Vio = 4 X19 — Xp — X%y5 — Xp. (13) 


Man veranschaulicht sich das Wirken des Operators am besten durch die 
Figuren 2a, b, c, welche die Koeffizienten enthalten, mit denen die jeweiligen 
Funktionswerte der Funktion x zu nehmen sind. . 
Fur die quadratische Form Q(x) ergibt sich der folgende Ausdruck: Man 
erganze zundchst die in den inneren Punkten gegebene Funktion x durch Hinzu- 
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fiigen von Randwerten, die alle Null sind, zu einer in allen Gitterpunkten er- 
klarten Funktion x. Ferner bezeichne man mit A, B die Endpunkte einer be- 
liebigen (im Inneren oder auf dem Rande gelegenen) Gitterstrecke. Dann ist 


Q(x) = DI (x4 — x5), (14) 


wobei tiber alle Gitterstrecken zu summieren ist. Diese Form ist positiv definit. 
Sie kann namlich nur verschwinden, wenn zwei durch eine Gitterstrecke ver- 
bundene Werte von x gleich sind. Da die Randwerte Null sind, ist dies nur 


ae eas eae oS 


2 4 { { 2 
e e . e e ° 
{ { 
-4 e 
Im Inneren 
i 4 Null { 
+4 ~4 
+h 1 +4 | 
-4 i} 
a b c Deere ene mt * 
Innere Punkte Punkt 10 Punkt 5 a Ses 
Fig. 2 Fig. 3 


moglich, wenn x identisch Null ist. (14) vertritt natiirlich das Dirichlet-Integral 
der kontinuierlichen Theorie. 

Der erste Schritt zur Lésung unseres Randwertproblems besteht nun in der 
Wahl einer Versuchsgitterfunktion v und der zugeh6rigen Restduenrechnung 
(Bestimmung der Residuengitterfunktion 7). Dabei schreibt man natiirlich die 
25 Gleichungen nicht hin, sondern bildet entsprechend der Formel (3’) die 
Funktion Dv unter Benutzung der Figuren 2a, b, c. Man kann diese Rechnung 
durch die folgende Summenprobe kontrollieren. Es sei e diejenige Gitterfunk- 
tion, die in allen inneren Gitterpunkten den Wert 1 hat. Dann ist unter Verwen- 
dung von (3’) die Summe aller Residuen 


(r, e) = (Du + 1, e) = (Dv, e) + (t, e) = , De) + (2), (15) 


wobei die letzte Umformung aus der Greenschen Formel (9) folgt. Hierin ist 
(, e) die Summe der Konstanten in unseren 25 Gleichungen, also nach (10), 
(11) die negative Summe aller gegebenen Randwerte (mit Ausschlu8 der vier 
Quadratecken). Die Werte von De sind in Figur 3 angegeben. Die Summe der 
Residuen hangt also nur ab von den Funktionswerten von v in denjenigen Git- 
terpunkten, die in Figur 3 eingetragen sind. 

Fiir das Folgende brauchen wir noch Gitterfunktionen, die nur in einem 
einzigen Gitterpunkt von Null verschieden sind. Es sei zum Beispiel é die 
Funktion, die im Gitterpunkt Nr. 9 den Wert +1, aber in allen anderen Gitter- 
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punkten den Wert Null hat. Ist dann x eine beliebige Funktion und y = Dx, 
so gilt unter Beriicksichtigung von (12) 


(Dég, x) = (€9, Dx) = (€, Y) = Vo = 4 %q — %q — Xg — XQ — Xg - 


Vergleicht man nun im ersten und letzten Term dieser Gleichung die Koeffi- 
zienten der x;, so ergibt sich, daB Deg genau durch die Figur 2a gegeben ist, 
also im Punkt 9 den Wert 4, in den Nachbarpunkten den Wert —1 und sonst 
iiberall den Wert Null hat. Analog ist De, die Figur 2b und De; die Figur 2c. 

Man kann diese Bemerkung dazu benutzen, um die neuen Residuen 7’ nach 
einer Funktionskorrektur Av gemaB8 (4’) auszuwerten. Die Funktionskorrektur 
bestehe zundchst darin, daB nur der Funktionswert im Punkte 9 um einen 
Betrag « geandert werde. Dann ist 


Av=ae,, 
also 
f=7 + DAt=7 po De,. 


Man erhalt also die neue Residuenverteilung, indem man zur alten das «-fache 
der Figur 2a addiert. Eine Anderung von v um +1 im Punkt Nr. 9 bewirkt also 
eine Residuenaénderung von +4 im Punkte Nr. 9 und Residuenadnderungen — 1 
in den vier Nachbarpunkten. Aus diesem Grunde nennt man die Figuren 
2a, b,c auch Relaxationsfiguren. Wird v in mehreren Gitterpunkten simultan 
korrigiert, so muB man eben die einzelnen 
Residuenanderungen superponieren. 

In Figur 4 ist der Start einer Relaxations~ 
rechnung illustriert. Links von jedem Git- 
terpunkt sind Funktionswerte oder Funk- 
tionsanderungen angeschrieben, rechts wird 
der jeweilige Stand der Residuen protokol- 
liert. Es wurde mit einer Versuchsfunktion v 
begonnen, die sich verniinftig an die Rand- 
werte — die ja vorgegeben sind — anschlieBt, 
und die Residuenrechnung ausgefiihrt. 
Dann erfolgte ein erster Relaxationsschritt. 
Um das Residuum im Punkte Nr. 9 herab- 
zudriicken, wurde dort die Funktion um 
+8 geadndert. Dies ergibt nach der eben 
hergeleiteten Regel die eingetragenen neuen 
Residuen im Punkt Nr. 9 und seinen vier Nachbarn. Man miiBte nun sukzessive 
mit derartigen Punktrelaxationen weiterfahren, wobei natiirlich jeder Gitter- 
punkt mehrmals erfaBt werden muB. Um die Rechnung nach mehreren Schritten 
zu kontrollieren, summiert man die bisherigen Punktionsanderungen auf, bildet 


Fig. 4 
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So eine verbesserte Versuchsfunktion und fithrt mit ihr die Residuenrechung 
wieder durch. Ein geiibter Rechner, der ein gewisses Gefiihl fiir den Verlauf 
einer Potentialfunktion hat, kann im vorliegenden Beispiel von 25 Unbekannten 
mit diesem Einzelschrittverfahren durchkommen!). Im allgemeinen muB aber 
die Methode noch bedeutend verfeinert werden. Auch muB untersucht werden, 
wann man aufhoren darf, das heiBt, wie weit die Residuen herabgedriickt werden 
mussen, um die Lésung w des Dirichlet-Problems mit einer gegebenen An- 
naherung zu erreichen. 

Es ist klar, daB auch andere Randwertaufgaben (zum Beispiel diejenigen 
der Elastizitatstheorie) auf diese Weise angegriffen werden kénnen. 


2. Das allgemeine Verfahren der Relaxationsrechnung 
Wir kehren zum allgemeinen Fall der Gleichungen (1) und zur Sprache der 


Vektoren zurtick und bilden zu einem gegebenen Versuchsvektor v die quadra- 
tische Funktion 


F(o) = = Qe) + (9). (16) 


Q(v) ist die in (7) definierte quadratische Form Q(v) = (v, Dv) und / der Vektor 
der Konstanten. Ausgeschrieben: 


1 
F(v) = = Dd tite + DIL 0;. (17) 
~ (1, (2) 
Differenzieren wir das nach v,, so ergibt sich 


OF 
dv =D, Giete ela. 
(h) 


Dies ist aber nach (3) nichts anderes als das 7-te Residuum 7;: 


OF 
= — 1 
e Ou; (18) 
Wir konnen auch schreiben 
r =igrad tf; (19) 


was eben eine abgekiirzte Schreibweise fiir die Aussage bedeuten soll, dab die 
t-te Komponente von 7 gleich der partiellen Ableitung von F nach v, ist. Hieraus 


folgt zundachst: 


1) Rinzelschrittverfahren bedeutet: Bei jedem Relaxationsschritt wird die Versuchstunktion 
nur in einem einzigen Punkt geandert. Vgl. dazu auch A. WiLLErRs, M ethoden der praktischen Ana- 
lysis, 2. Aufl. (W. de Gruyter, Berlin 1950), S. 311. 
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Satz 1: Die Auflisung der gegebenen Gleichungen (1) ist vollstandig aqui- 
valent mit der Bestimmung des Minimums der quadratischen Funktion F). 

Denn um die Gleichungen aufzulésen, miissen die Residuen zu Null gemacht 
werden, und nach (18) bedeutet dies das Aufsuchen derjenigen Werte v,, fiir 
welche alle partiellen Ableitungen von F verschwinden, das heibt, fiir welche & 
stationar ist. Da die quadratische Form Q(x) als positiv definit vorausgesetzt 
ist, hat F genau eine stationare Stelle, und zwar ein Minimum. 

Satz 1 kann zur Begriindung des folgenden Relaxationsverfahrens dienen. 
Ausgehend von einem Versuchsvektor v = (vj, Up, ..., Up) und nach Wahl eines 
vom Nullvektor verschiedenen Gewichtsvektors p = (f,, po, ---, Py) Mache man 
die Korrekturen 


Av;=4; (20) 


der Versuchswerte, wobei A ein noch zu bestimmender Parameter ist. Schreibt 
man (20) vektoriell: 


Av=iA?, (21) 
so betrachtet man also die lineare Schar von Versuchsvektoren 
v=v+Av=v+Ap. (22) 


Bildet man mit jedem v' die quadratische Funktion F, so wird F eine Funktion 
F(A) von A. Man bestimme nun A so, daB F (als Funktion der einzigen Variablen 
A betrachtet) minimal wird. Es handelt sich hier im Grunde genommen um eine 
Anwendung des Ritzschen Gedankens, ein Minimum angenahert dadurch zu 
bestimmen, daB man die Konkurrenz auf eine leicht zu handhabende lineare 
Schar beschrankt. Es ergibt sich so unter Benutzung von (18) 


dF (A) OF dv; : : 
as ct Ou, dA D> hh (Deo: (23) 
\ : (?) 


Hierin ist 7’ der Residuenvektor, der zum Versuchsvektor v’ gehort. Mit anderen 
Worten: . 

Satz 2: Nach Durchfiihrung eines Relaxationsschrittes muB der neue Resi- 
duenvektor orthogonal auf dem Gewtchtsvektor stehen. 


Jetzt muB nur noch die Orthogonalitaétsbedingung ausgewertet werden. 
Setzt man (4’) in (23) ein, so folgt 


(r + DAv, p) = (r +A Dp, p) = (r, p) +A (Dp, ) = 0, 


1) Bekanntlich kann man viele Randwertaufgaben der mathematischen Physik direkt als 
Variationsprobleme formulieren (Energiemethode), so da8 sich beim Ubersetzen in die Differenger 
rechnung direkt die Funktion *(v) ergibt, welche extremal gemacht werden muB. Die zugehorigen 
linearen Gleichungen ergeben sich dann nach (18) durch Differenzieren von F. Dieser Weg ist zu 
empfehlen, da infolge der physikalischen Bedeutung von F (Energie) meistens von vorneherein fest- 
steht, daB der homogene Teil, das heiBt die quadratische Form Q, positiv definit ist. 
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also 


pei wale. Bhat” 

(DP, p) ee 
damit ergeben sich folgende Regeln fiir einen allgemeinen Relaxationsschritt: 
E's set v ein Versuchsvektor, r der zugehirige Residuenvektor. Nach Wahl eines 
Gewichtsvektors p mache man die Korrektur 


AveSieg Maat) 7D. (25) 


Die neuen Residuen sind dann gegeben durch 
Yost A Dp, (26) 


Aus der Herleitung folgt, daB der Wert der quadratischen Funktion F bei 
jedem Relaxationsschritt sicher abnimmt. Denn fiir die zweite Ableitung von 
F(A) ergibt sich aus (23) und (26) 


ae = (Go?) = @ DP) = 0) > 0. (27) 


Diese Ableitung ist also eine positive Konstante und somit die graphische Dar- 
stellung von F(A) eine nach oben gedffnete Parabel. 

Bei jedem Relaxationsschritt nahert man sich also dem schlieBlich gesuch- 
ten Minimum von /’(v). Hierin liegt der Beweis fiir die Konvergenz dieses Re- 
laxationsverfahrens gegen die Lésung des gegebenen 
Gleichungssystems?). Fiir die systematisch gelenkten 
Relaxationsverfahren (wie das Einzelschrittverfah- 
ren oder das Verfahren des starksten Abstiegs) kann 
diese Bemerkung zu einem mathematisch exakten 
Beweis ausgestaltet werden dafiir, daB die Versuchs- 
vektoren gegen den Lésungsvektor konvergieren. 
Die auseinandergesetzten Verhaltnisse lassen sich 
im Falle von » = 2 Gleichungen geometrisch illu- 
strieren (Figur 5). Im Koordinatensystem 1,, v, ist 
die quadratische Funktion F(v) durch ihre Niveau- 
linien dargestellt. Es sind dies konzentrische Ellip- i 


sen. Das gemeinsame Zentrum ist das Minimum, 
also sind seine Koordinaten 2, vu, die Lésungen des 
Gleichungssystems. Es wurde ein Versuchspunkt v 
eingezeichnet; der zugehdrige Residuenvektor 7 steht nach (19) senkrecht auf 
der Niveaulinie. Ein Relaxationsschritt besteht dann darin, daB nach Wahl 


1) R. von Mises und H. PoLLaczeK-GEIRINGER, Z. angew. Math. Mech. 9, 58-77, 152-164 
(1929); vgl. auch L. Cottarz, Zur Herlettung von Konvergenzkriterien fiir Iterationsverfahren bei 
linearen Gleichungssystemen, Z. angew. Math. Mech. 30, 278-280 (1950) ; vgl. ferner G. TreMPLE, dives 
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eines Richtungsvektors # von v aus in der Richtung p vorwartsgegangen 
wird bis zu dem Punkt v’, in dem F auf diesem Wege minimal wird. Dieser 
Punkt befindet sich offenbar dort, wo der Weg eine Niveaulinie beriihrt. Dort 
steht der neue Residuenvektor als Gradient senkrecht auf der Niveaulinie, also 
senkrecht auf p (Satz 2). 

Beim Einzelschrittverfahren wird jeweilen nur v, oder vy geandert, das heiBt, 
man benutzt nur zwei Relaxationsrichtungen # wahrend des ganzen Relaxa- 
tionsprozesses, namlich die Parallelen zu den Koordinatenachsen. 

Die Sache wird noch iibersichtlicher, wenn man die Figur 5 so in eine Figur 6 
affin abbildet, daB die Ellipsen zu Kreisen werden. Ein Relaxationsschritt be- 
deutet dann einfach, daB man den Naherungspunkt wv 
orthogonal auf eine gewahlte Gerade « projiziert, die 
durch das Zentrum lauft. Wahlt man nur zwei feste 
Gerade « und f, auf die man abwechselnd projiziert — 
dies ist zum Beispiel beim Einzelschrittverfahren dér 
Fall —, so verlauft die Relaxation auf dem Zickzackweg 
der Figur 5. Wir bezeichnen den kleineren Winkelraum 
zwischen « und f als Kdfig, weil man leicht feststellt, 
daB die Relaxation — gleichgiiltig, wo sie beginnt — ein- 
mal in diesen Kafig eintritt und alsdann dauernd in ihm 

Fig. 6 bleibt. Je schmaler dieser Winkelraum ist, um so lang- 

samer konvergiert das Verfahren gegen den Lésungs- 

punkt. Das Auftreten von Kafigen ist eine allgemeine Erscheinung bei Rela- 

xationsverfahren und sehr unerwiinscht. Es bewirkt, daB die Relaxation am 
Anfang flott vorwartsgeht, aber dann immer weniger ausgiebig wird. 

Im Falle von m > 2 Gleichungen kann man die Verhaltnisse natiirlich analog 
in einem n-dimensionalen Raum veranschaulichen. 

Die dargebotenen Uberlegungen sollen nun noch auf das Dirichlet-Problem 
angewendet werden. Die in Figur 4 gezeigte Punktrelaxation (Einzelschritt) er- 
scheint nun wie folgt als Spezialfall der allgemeinen Theorie: Da die Versuchs- 
funktion nur im Punkte 9 geindert werden soll, ist die Gewichtsfunktion p 
gleich ég, also 1 im Punkte Nr. 9 und sonst iiberall Null. Nach Satz 2 muB die 
Relaxation so durchgefiihrt werden, daB die neue Residuenfunktion 7’ ortho- 
gonal zu éy ist, das heiB®t, daB der neue Residuenwert im Punkte 9 Null wird. Es 
ergibt sich also die 

Regel: Bet einer Punktrelaxation soll das Residuum im betreffenden Punkt 
vollstindig beseitigt werden. Entsprechend der Relaxationsfigur 2a verteilt es sich 
dabei gleichmapig auf seine vier Nachbarpunkte?) (vergleiche Punkt 2 in Figur 4). 

Es ist leicht festzustellen, wann man bei diesem Einzelschrittverfahren des 
Dirichlet-Problems in einem Kafig sitzt. Dies ist dann der Fall, wenn die Resi- 


1) Die analoge Regel gilt natiirlich auch fiir das Einzelschrittverfahren bei beliebigen Gleichungs- 
systemen. 
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duen in den inneren Punkten alle dieselben Vorzeichen haben. Nehmen wir 
etwa an, sie seien alle = 0. Da sich bei jeder Punktrelaxation das Residuum 
im betreffenden Punkt auf die vier Nachbarn gleichmaBig verteilt, sind auch 
nachher wieder alle Residuen = 0, das hei®t, die Relaxation verlaiuft weiterhin 
dauernd im Kafig. Der Leser mag sich an einem einfachen numerischen Beispiel 
tiberzeugen, daB dies auch hier eine starke Verminderung der Konvergenzge- 
schwindigkeit bedeutet, so daB der positive Residuenberg mit dem Léffel statt 
mit einer Baggermaschine abgetragen wird. Andererseits besteht die Gefahr, 
da8 wie in Figur 6 der RelaxationsprozeB — gleichgiiltig, mit welcher Residuen- 
verteilung er startet — einmal in den Kafig eintritt oder sich wenigstens groBere 
Teilgebiete positiver Residuen (Teilkdfige) bilden, zum Beispiel in der linken 
oberen Ecke und am oberen Rand der Figur 4. Es bestehen folgende Méglich- 
keiten, um einen derartigen Teilkafig abzubauen und damit die unbedingt not- 
wendige Konvergenzbeschleunigung zu erreichen. 


a) Blockrelaxation 


In Figur 7 ist ein Teil des quadratischen Gitters im Grundgebiet gezeichnet, 
der einen Kafig bilden soll, und zwar seien alle Residuen in den schwarzen Git- 
terpunkten positiv. Dann wahlen wir als Gewichtsfunktion # fiir einen allge- 
meinen Relaxationsschritt (25) die Funktion, welche in allen schwarzen Punk- 
ten gleich 1, auBerhalb aber Null ist. In Figur 7 wurden die Werte der Funktion 


Dp angeschrieben und (f, Dp) = Q(p) = 22 
ausgerechnet. [Wegen (14) ist Q(f) gleich der 
Anzahl der gestrichelten Strecken.] Die ganze 
Figur 7 bezeichnen wir als «Block». 

Nach der Regel (25) ist nun ein solcher 
Block folgendermaBen anzuwenden: Man er- 
halt A, indem man die negativ genommene 
Summe der Residuen in den schwarzen Punk- 
ten durch 22 dividiert. Die Versuchsfunktion 
wird in allen schwarzen Punkten um denselben 
Wert 4 korrigiert. Die dadurch bewirkten 
Residuenanderungen ergeben sich, indem man 
die Zahlen der Figur 7 mit A multipliziert. 

Nach Satz 2 hat dieser Relaxationsschritt 


x ap eh 
e 9 
i] 
1 


ie. 


0 0 


folgende Wirkung: Die Summe der neuen Residuen in den schwarzen Punkten 
wird Null. Da aber nach Figur 7 die Residuen in den inneren Punkten nicht 
gedndert werden, bedeutet dies, daB auf dem dick gezogenen Rand der Figur 7 
positive und negative Residuen hergestellt werden, im Innern aber ein posi- 
tiver Kafig bleibt. Dieser mu8 dann durch einen kleineren Block angegriffen 


werden. 
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b) Schetbenrelaxation 


Man kann die Blockmethode verallgemeinern, indem man fiir die Gewichts- 
funktion p eine beliebige, in den schwarzen Punkten subharmonische Funktion 
nimmt, die auBerhalb Null ist. (Subharmonisch bedeutet Dp > 0.) Nach (26) 
aindern sich namlich dann alle Residuen bei diesem Relaxationsschritt im glei- 
chen Sinn. Eine derartige Gewichtsfunktion nennen wir eine Schezbe. Unter 
den vielen méglichen Methoden, Scheiben zu bauen, seien die beiden folgenden 
erwahnt: 

1. Das Gebiet der schwarzen Gitterpunkte, fiir welche eine Scheibe kon- 
struiert werden soll, sei ein Rechteck. Dann lése man das Eigenwertproblem 


Dp = up, in jedem schwarzen Punkt, (28) 
~=0 auBerhalb (29) 


(schwingende Membran) und nehme fiir # die Grundschwingung; diese hat be- 
kanntlich keine Knoten, also in allen schwarzen Punkten dasselbe Vorzeichen, 
somit ist nach (28) tatsachlich # subharmonisch. 

Dieses Eigenwertproblem ist durch Fourier-Entwicklung leicht zu lésen. 
Es ergibt sich zum Beispiel fiir ein Quadrat von 5 - 5 Gitterpunkten im Koor- 
dinatensystem der Figur 8 


p(x, y) =sin =*'sin =>, (30) 
und der Eigenwert dieser Grundschwingung ist 
Je oat 
\ = 8sin® Rig: 0,536 . 


In Figur 8 sind die Werte von # mit drei Stellen 
83 tabelliert. Noch ein Beispiel zur Anwendung der 
50 Scheibe: die gegebene Residuenverteilung bilde 
einen Teilkafig, der sich ungefahr auf ein quadra- 
tisches Gebiet von 5-5 Gitterpunkten erstreckt, 
und zwar so, daB das Residuum in jedem dieser 
Punkte ungefahr + 100 betragt. Die Anwendung 
der Scheibe ergibt dann nach Regel (25) im 
(b, Dp) = (P,P) = 90000 uw = 48240 Mittelpunkt des Quadrats eine Funktionskorrek- 
Fig. 8 tur von etwa —(289). Demgegeniiber liefert eine 
einzelne Punktrelaxation in diesem Mittelpunkt 
nur eine Funktionskorrektur von —100 :4 = —25 zur Beseitigung des Residu- 
ums --100, also einen mehr als zehnmal zu kleinen Wert! 
2. Kleinere Scheiben, die nur wenig Gitterpunkte enthalten, kénnen direkt 
durch Relaxieren konstruiert werden. Man geht aus von einer Versuchsfunktion 


Sra 
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fir # und bildet Dp. Durch systematisches Korrigieren der Versuchsfunktion 
in Einzelschritten beseitigt man zundchst alle negativen Werte von Dp und 
gleicht dann die positiven Werte von Dé aus, so daB zum SchluB Dp iiberall 
fast iibereinstimmende positive Werte hat. 

Gr6Bere Scheiben bildet man, indem man innen in eine kleinere bekannte 
Scheibe ein Feld von Versuchswerten einsetzt und langs der entstandenen Naht 
durch Relaxation ausgleicht. 

Scheiben sind zum Aufbrechen von Kafigen im allgemeinen wirksamer als 
Blécke, da sie im Gegensatz zu diesen auch das Innere des Kafigs angreifen. 
Es ist sogar gelegentlich anzuraten, Teilkafige absichtlich herzustellen, um 
Scheiben anwenden zu kénnen. Dies kann zum Beispiel so geschehen, da in 
einem Teilgebiet langere Zeit — etwa nach dem Gesamtschrittverfahren von 
Abschnitt 4 — relaxiert wird. Wir werden spater zeigen, daB dann die Residuen- 
verteilung asymptotisch die Form der oben eingefiithrten Grundschwingung 
dieses Teilgebietes bekommt, so daB eine Scheibe vom Typus 1 wunderbar 
«paBt». 

Es ist zweckmabig, fiir einen gegebenen Operator D eine Sammlung von 
Blocken und Scheiben anzulegen. Zum Beispiel kann eine solche Sammlung 
von Blécken und Scheiben fiir den Operator der Figur 2a dann zur Lésung jedes 
Potentialproblems Verwendung finden, gleichgiiltig, wie das Grundgebiet aus- 
sieht oder wie die Randbedingungen lauten. 


3. Simultane Relaxation 


Zuriickkehrend zum Beginn des Abschnitt 2 k6nnen wir den dortigen Ansatz 
(21) wie folgt verallgemeinern: An Stelle des einzigen Gewichtsvektors p wahlen 
wir jetzt m Gewichtsvektoren fy, 9, .--, Pm, die linear unabhangig sein sollen’). 
Dann machen wir fiir den Korrekturvektor des Versuchsvektors v den Ansatz 


Ay = Ay Py de Pa + + Am Pm (31) 
mit den jetzt zu bestimmenden Parametern Aj, dy, ..., A. Wir betrachten also 
die lineare Schar von Versuchsvektoren 

v=vt+Av=o+ DA pe. (32) 
kel 


Wie in (23) ergibt sich fiir die Ableitung der zu minimalisierenden Funktion 
F(v') nach A;: 
Se Spe 33 


wobei 7’ der zu v’ gehérige Residuenvektor ist, also 


1) Indizes bedeuten von jetzt an Nummern von Vektoren und nicht mehr wie friher die Kom- 
ponenten eines einzelnen Vektors. 
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Satz 3: Bei simultaner Relaxation mit mehreren Gewichten muf der neue 
Residuenvektor orthogonal zu saémtlichen Gewichtsvektoren sein. 
Die Auswertung dieser Orthogonalitaétsbedingung ergibt nach (4’) und (31) 


vy =r+ DAv=rt+ DAEOP ES 
bas 


(1', Ps) = (1, Bi) + 2 An (bi, Db,) = 0. 


Man erhalt also zur Bestimmung von Aj, Ag, ..., A, das lineare Gleichungs- 
system 


m 


dD. (bir Dix) An + (1, Bi) = 0. (34) 


k=1 


Dabei ist 7 der Residuenvektor, der zum Versuchsvektor v gehdrt. Wegen der 
Greenschen Formel (9) ist dieses System symmetrisch, fallt also seinerseits wie 
das Ausgangssystem (1) unter die Theorie von Abschnitt 1 und kénnte daher 
durch Relaxationsrechnung aufgelést werden. Die zu ihm gehirige quadratische 
Form ist nach (7) 


>? (Di, Dp) Ai Ax a pa A; P: , D Dm ps) ’ 


(4, k) k=1 


und dies ist nichts anderes als unsere alte Form 


OA, Pr at Ay pe ret Aw Pm) : 


Sie ist daher immer gréBer als oder gleich Null und kann nur verschwinden, wenn 


Ay by + Ag Po + ess An Dim =0, 


was wegen der Unabhiangigkeit der ; nur méglich ist, wenn 
Ay =A, =A, =A, ge Am =(. 


Damit ist gezeigt, daB das Gleichungssystem (34) die in Abschnitt 1 gestellte 
Definitheitsforderung erfiillt. Sind die Gleichungen (34) aufgelést, so findet man 
die Funktionskorrekturen und neuen Residuen aus 


nt 


Av= Saudr, Y=rt D> 4x DP, . (35) 
k=1 k=1 


Im allgemeinen wird man natiirlich die Zahl der Gewichte nicht zu groB und 
speziell kleiner als die Zahl 1 der Unbekannten in (1) wahlen. Die Methode der 
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simultanen Relaxation besteht also — grob gesprochen — darin, daB das Glei- 
chungssystem (1) mit Hilfe von Gleichungssystemen mit weniger Unbekannten 
gelost wird. Wir geben im folgenden drei Anwendungen der Methode. 


a) Verbesserung der Blocktechnik 


Wir haben in Abschnitt 2 festgestellt, daB ein einzelner Block p (Figur 7) 
gewohnlich nicht gentigt, um einen Kafig abzubauen, sondern daB anschlieBend 
noch mehrere Blécke im Innern des ersten anzuwenden sind. Es soll etwa in 
Figur 7 auBer # = fp, noch ein weiterer Block #, auf das Quadrat ABCD gelegt 
werden, so daB also die Gewichtsfunktion £, in jedem der neun Punkte dieses 
Quadrats den Wert +1, sonst aber iiberall den Wert Null hat. Statt nach- 
einander sollen die Blécke nun simultan angewandt werden, was viel vorteil- 
hafter ist. Man berechnet sofort 


(1, Dp;) = 22, (bz, DPy) = 1, (Pe, De) = 12, 
so daB die Gleichungen (34) lauten: 
2244+ 4,+ (7,2:))=9, A44+124,+ (7, b) =0. 


Hierin ist (vr, #,) die Summe der Residuen in den schwarzen Punkten der 
Figur 7 und (7, #,) die Summe der Residuen in den neun Punkten des Qua- 
drates ABCD. Die Wirkung dieser simultanen Relaxation ist nach Satz 3 die, 
daB sowohl am Rande des Quadrats A BCD als auch in den schwarzen Punkten 
auBerhalb dieses Quadrates Residuen beiderlei Vorzeichens entstehen. Hatte 
man diesmal den Block f, ganz ins Innere von 4, gelegt, so ware (p,, Dp,) = 0, 
das heiBt, es kame auf dasselbe heraus, wenn die beiden Blécke simultan oder 
nacheinander angewendet wiirden. 


b) Konjugierte Gewichte 


Wir nennen zwei Gewichtsvektoren #, und f, konjugiert, wenn (p, Dp,) = 0 
ist. Wegen der Greenschen Formel ist dann auch (,, Dp,) = 0, so daB also diese 
Relation zwischen zwei Gewichten symmetrisch ist. Wir haben eben einen 
solchen Fall angetroffen. Die simultane Relaxation mit Gewichtsvektoren 
py Po +++» Pm die paarweise zueinander konjugiert sind’), ist besonders 
einfach, indem man an Stelle der Gleichungen (34) fiir Ay, Ag, -.., Am die expli- 
ziten Formeln hat: ag 

Vy Pi 
A ao (Pi DP;) oe 


1) Sind alle p; von Null verschieden, so sind sie von selbst linear unabhangig. 


16 EDUARD STIEFEL ZAMP 


Gelingt es fiir ein gegebenes Problem (1) mit » Unbekannten auch 7 von Null 
verschiedene und paarweise konjugierte Gewichte zu finden, so fithrt eine ein- 
zige simultane Relaxation mit diesen Gewichten zur exakten Lésung des Pro- 
blems. Denn der neue Residuenvektor nach dieser Relaxation, 7’, steht nach 
Satz 3 orthogonal auf jedem der x linear unabhangigen Gewichtsvektoren, ist 
also der Nullvektor. 

Beim Dirichlet-Problem fiir ein gegebenes Grundgebiet kann man in einem 
rechteckigen Teilgebiet konjugierte Gewichtsfunktionen wie folgt konstruieren: 
Man lise wie bei der Scheibenkonstruktion das Eigenwertproblem (28) 


Dp=up in den Gitterpunkten des Teilgebietes, 


37 
p=0 auBerhalb des Teilgebietes 7 


durch Fourier-Entwicklung nach den beiden Koordinatenrichtungen der Figur 8 
und nehme als Gewichte #,, fo, ... die sukzessiven Eigenfunktionen, geordnet 
nach wachsenden Eigenwerten. Diese Funktionen sind orthogonal, also wegen 
(37) auch konjugiert. Die Berechnung der Koeffizienten (36) fiir die Simultan- 
relaxation ist dann nichts anderes als die zweidimensionale Fourier-Analyse der 
gegebenen Residuenverteilung im vorliegenden Rechteck, die Berechnung der 
Korrekturen (35) eine darauffolgende Synthese. Man kann also die sehr weit 
ausgebildeten Rechenverfahren der trigonometrischen Interpolation benutzen. 
Die Wirkung der Simultanrelaxation besteht dann darin, daB , als Scheibe 
einen etwa vorhandenen Kafig beseitigt, wahrend die héheren Eigenfunktionen 
dann noch vorhandene Rauhigkeiten glatten. Fiir ein Quadrat aus vier Gitter- 
punkten erhalt man so die in Figur 9 gezeichneten vier konjugierten Gewichts- 


imal io eh pat 
/ 1 0 ~4 4 0 4 ~4 
Fig. 9 


funktionen. Eine simultane Relaxation mit diesen vier Gewichten macht die 
Residuen in den vier Quadratecken exakt zu Null. 


c) Ausrelaxieren eines Teilgebietes 


Es sei etwa beim Dirichlet-Problem wieder ein Teilgebiet von m Gitter- 
punkten gegeben. Denkt man sich fiir einen Moment m Gewichtsfunktionen 
gewahlt, von denen jede in genau einem Gitterpunkt des Teilgebiets den Wert 1 
hat, sonst aber tiberall den Wert Null. Dann besagt Satz 3, daB die zugehérige 


Simultanrelaxation die Residuen in simtlichen m Gitterpunkten des Teilgebiets 
zu Null machen muB. 
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Ohne die allgemeine Theorie zu benutzen, kann man die Gleichungen fiir 
die Funktionskorrekturen im Teilgebiet direkt hinschreiben. Ist zum Beispiel 
in Figur 1 das Teilgebiet die Strecke der Punkte Nr. 6, 7, 8, 9, so lauten diese 
Gleichungen: 

4 Avg— Av, +7%=0, 
— Ave +4Av,— Avg +7, =0, | 
= As, 4Ay,— Av 7,=0, + ce 
— Avet+4Avy+7,=0. | 


Hierauf kann man folgende Relaxationstechnik griinden: Man wahlt eine ein- 
fache Punktfigur (zum Beispiel ein Quadrat aus 3 - 3 Gitterpunkten), bringt 
diese Figur an verschiedene Stellen im Grundgebiet und beseitigt jeweilen die 
Residuen im Inneren der Figur durch Auflésen der Gleichungen vom Typ (38). 
Das Einzelschrittverfahren (Punktrelaxation) ist natiirlich derjenige Spezial- 
fall, wo die Figur aus einem einzigen Punkt besteht. 

Das haufige Auflésen der Gleichungen wird man sich erleichtern, indem man 
ein fiir alle Male die inverse Matrix der Gleichungsmatrix bildet. Wahlt man 
Punktfiguren hoher Symmetrie, so bietet die Theorie der Charaktere der zuge- 
hérigen Symmetriegruppe Méglichkeiten, um die Gleichungen in Systeme mit 
weniger Unbekannten zu zerfallen. 

Im Institut fiir angewandte Mathematik der ETH. in Ziirich wurden die 
elastischen Spannungen in einer Gewichtsstaumauer unter Beriicksichtigung 
der elastischen Einbettung in den Baugrund so bestimmt?). [Aufldsung der 
biharmonischen Gleichung fiir die Airysche Spannungsfunktion durch Rela- 
xation, Punktfigur = Quadrat von 4 - 4 Punkten, Auflésung der 16 Gleichungen 
vom Typus (38) vollautomatisch auf einer programmgesteuerten Rechenma- 
schine. | 


4. Methode des starksten Abstiegs 
a) Gesamtschritiverfahren 


Wahrend die bisherigen Methoden der Intuition des Rechners betreffend 
der Wahl des nachsten Relaxationsschrittes weitgehenden Spielraum lassen, 
kommen wir jetzt zu den Verfahren, die zwangslaufig fortschreiten und sich 
daher fiir programmgesteuerte Rechenautomaten eignen. Diese Verfahren ge- 
hoéren alle in den Abschnitt 2, arbeiten also pro Schritt mit einem einzigen 
Gewichtsvektcr # und nicht simultan mit mehreren Gewichten. Nach jedem 
Schritt wird der nichste Gewichtsvektor vorgeschrieben und ist also nicht wie 
bisher noch wahlbar. 

1) Ein Bericht daritiber erscheint demnachst in den Comptes rendus des Colloques internationaux 


du Comité national de recherches scientifiques, Paris, Colloque de l’Institut Blaise-Pascal, 8. bis 
13. Januar 1951. 


ZAMP III/2 
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Wir werfen einen Blick auf die Figur 5. Ausgehend von einem Versuchs- 
punkt v sollen wir in der Richtung des Gewichtsvektors p vorwartsgehen, so — 


daB unsere quadratische Funktion / abnimmt. Nun nimmt aber /’ (wenigstens 


in der Umgebung von v) dann am stirksten ab, wenn wir senkrecht zu den — 


Niveaulinien, also in der entgegengesetzten Richtung zu grad F, fortschreiten. 
Ist y der zu v gehérige Residuenvektor, so bedeutet dies wegen (19) die Wahl 


poor, (39) 


Das Verfahren des stiirksten Abstiegs besteht also darin, daB der neue Gewichts- 


vektor entgegengesetzt gleich dem Residuenvektor gewahlt wird. Satz 2 sagt — 


dann: 

Satz 4: Beim Verfahren des stirksten Abstiegs steht nach Durchfiihrung eines 
Relaxationsschrittes der neue Residuenvektor orthogonal auf dem alten. 

Die Formeln (25) ergeben die Regeln: 

Es sei v ein Versuchsvektor, 7 der zugehérige Residuenvektor. Dann betragt 
die Korrektur 


(os =A eeniteeen as (40) 


und die neuen Residuen: 
Y =r—A Dr. (41) 


Der Quotient (40) spielt eine groBe Rolle beim Problem, die Eigenwerte der in 
Abschnitt 1, Formel (7), eingefiihrten quadratischen Form Q(x) = (x, Dx) zu 
bestimmen. Man bezeichnet gewohnlich 


R(x) = 24) (42) 
als Rayleighschen Quotienten des Vektors x. Sein Wert liegt immer zwischen 
dem gr6Bten und dem kleinsten Eigenwert der Form Q(x). Wir stellen also fest, 
daB unser Korrekturfaktor 4 der reziproke Rayleighsche Quotient des Residuen- 
vektors ist. 

Haben wir zum Beispiel beim Dirichlet-Problem einen Kafig iiber das ganze 
Grundgebiet hinweg, also tiberall positive Residuen, so liegt diese Residuen- 
verteilung in der Nahe der Grundschwingung der Membran, das heiBt in der 
Nahe der ersten Eigenfunktion des Eigenwertproblems Dx = x. Der Ray- 
leighsche Quotient von y ist somit in der Nahe des kleinsten Eigenwertes [ty 
und daher 4 grof, das heiBt nahe am maximalen Wert 1/1, den A iitberhaupt an- 
nehmen kann. Die Funktionskorrekturen sind nach (40) auch groB, was ja zur 
Beseitigung des Kafigs nétig ist. Haben hingegen die Residuen stark abwech- 


ye ; “A a7 . es . o . 
selnde Vorzeichen, so sind wir cher in der Nahe einer héheren Eigenfunktion 
und 4 wird klein. 
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Unter Gesamtschrittverfahren') verstehen wir folgende Variante der Methode 
des starksten Abstiegs: Man rechnet bei diesem Verfahren wie bei (40), (41) 
nach den Formeln 


Av=-—Ar, vr =r—ADr, (43) 


nimmt aber fiir 2 nicht den reziproken Wert des Rayleighschen Quotienten, 
sondern eine ein fiir allemal fest gewdhite Konstante, die bei allen Relaxations- 
schritten verwendet wird. Man erspart sich damit die Berechnung der skalaren 
Produkte (Dr, 7) und (r, 7) bei jedem einzelnen Schritt. Da es sich hier um eine 
Vorschrift handelt, die von unserem Grundprinzip des Minimalisierens der 
quadratischen Funktion / abweicht, miissen wir neu untersuchen, ob das Ver- 
_ fahren gegen die Lésung des Problems konvergiert. Wir brauchen zu diesem 
Zweck die Eigenwerte 


My Sfp Shs Ss Shy 
und die zugehdérigen Eigenvektoren 


Riven ont 


des Eigenwertproblems Dx = wu x oder — was auf dasselbe hinauslauft — der 
quadratischen Form Q(x). Wir starten die ganze Relaxationsrechnung mit 
einem Versuchsvektor vy und entwickeln den zugehérigen Residuenvektor 7, 
nach Eigenvektoren: 


To = Oy Hy + Oy Hy ts + Oy Xp. (44) 


Dann ist der Residuenvektor nach k-maliger Anwendung des Gesamtschritt- 
verfahrens: 


V5 Oy (1 — A pay)® y+ Oty (1 — A pe)* ag + oe (1 —A a)? %,, (45) 


was man leicht durch vollstandige Induktion wie folgt verifiziert: Aus (45) 
folgt 
Te pa — Ve = — hy (LA y)PA oy yo Oy (LA fin)® A en Xn» 


und dies ist wegen Dx, = mu; x; gleich 


ee = A (LA) * Daley — o> — A(L—A fen)*® Doty %_ = —A Dry, 
1) In leichter Verallgemeinerung einer in der Literatur eingeburgerten Terminologie (vgl. Wi1L- 
LERS, 1. c., Seite 304). Das Verfahren samt den zugehorigen Konvergenzuntersuchungen wurde wohl 
zum ersten Male ausfiihrlich geschildert von L. F. Ricuarpson, The Approximate Arithmetical 
Solution by Finite Differences of Physical Problems..., Phil. Trans. Roy. Soc. London [A], 210, 


307-357 (1911). 
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also 
Peery == 1, —ADr, 


in Ubereinstimmung mit der Vorschrift (43) fiir jeden Relaxationsschritt. Damit 
nun das Verfahren konvergiert, also die Residuen gegen Null streben, muB jede 
der Klammern in, (45) dem Betrag nach kleiner als 1 sein. Dies ergibt als not- 
wendige und hinreichende Konvergenzbedingung 


0<A< - : (46) 
Je nach Wahl von / zeigt das Verfahren verschiedenes Konvergenzverhalten. | 
Es soll dies noch etwas diskutiert werden, wobei zweckmaBig 4 = 1/A eingefiihrt — 
wird, so daB (45) lautet 


j HM, \k Ha\ki es My \k 
l—— Hy (1 = | Xy + Xo (1 = £2) NXg ry ia a che + An (1 == i Xn - (47) | 

Fall 1: u, <u. Dann strebt jeder Koeffizient in (47) monoton gegen Null, — 
und es ergibt sich asymptotisch | 


Vp Oy (1 pO Y Xj (48) | 


(falls nicht zufallig der erste Eigenwert jz, mehrfach ist oder der Fourier-Koeffi- 
zient «, verschwindet). Ist speziell 4 = «,, so wird der Beitrag der héchsten 
Eigenfunktion an die Residuenverteilung sofort beseitigt. 

Fall 2: (uy + Un)/2 << fu, . Die asymptotische Formel (48) gilt noch; — 
die Koeffizienten in (47) mit “, << ~ nehmen monoton und diejenigen mit 
fu < py oszillierend ab. 

Fall 3: w = (My + Mp)/2. Hier gilt statt (48) 


Tem (1 7, . [oy % + (—1)* a, Xn] . (49) 


bl 


Fall 4: u,/2 < uw < (M+ up)/2. Jetzt ist | 


re ~ tty (1 Hn. \* sti. | (50) 


Aus dieser Diskussion ergibt sich deutlich, daB wir zwei Typen von Relaxations- 
problemen unterscheiden miissen. 

Typus I: Probleme mit y4,,/2 S jy. In diesen Fallen hat man es in der Hand, 
durch geeignete Wahl von yw den EinfluB jeder Eigenfunktion rasch zu elimi- 
nieren. (Man wahle eben « gleich dem betreffenden Eigenwert.) Fiir die prak- 
tische Durchfiihrung wird man also jz das Intervall der Eigenwerte uw, < bsp 
durchlaufen lassen. 7 el 
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L'ypus II: Probleme mit uw, S m,/2. Hier bleibt der Einflu8 der ersten 
Eigenfunktion bis zum SchluB merklich, und in den Fallen 1, 2,3 wird die 
Residuenverteilung sogar immer mehr dieser Eigenfunktion ahnlich). Wir 
finden also hier die friiher an speziellen Beispielen besprochene Erscheinung 
des Kdafigs wieder und werden definieren: Unter einem Kdafig versteht man eine 
Residuenverteilung, welche ungefihr die Gestalt der ersten Eigenfunktion hat. Wir 
erkennen, daB die Residuenverteilung im Verlauf der Relaxation die Tendenz 
hat, diese Kafigform anzunehmen und dafB der Abbau des Kafigs nur langsam 
erfolgt. Neben der ersten Eigenfunktion sind natiirlich auch die Beitrage der- 
jenigen Eigenfunktionen gefahrlich, deren Eigenwerte unterhalb ,,/2 liegen. 
Bei den Problemen vom Typus II wird man mit mw hart an die Konvergenz- 
grenze 1/2 herangehen, also Fall 4 wahlen. Um Kafige aufzubrechen, empfiehlt 
es sich sogar gelegentlich, fiir einige Relaxationsschritte  < m,/2 zu wahlen 
(Methode der Uberkorrektion). Dies bewirkt allerdings, daB die Beitrage der 
hoheren Eigenfunktionen wieder aufgeschaukelt werden?). Beim Verfahren des 
starksten Abstiegs wird A als Rayleighscher Quotient automatisch reguliert. 

Was liefern diese Erkenntnisse fiir das Dirichlet-Problem in einem beliebigen 
Grundgebiet, das quadratisch unterteilt wurde? Ist wu; ein Eigenwert dieses 
Problems, so ist auch (8 — w,) ein solcher, man kann also die Eigenwerte 
paaren. Um dieses einzusehen, nehmen wir die «Schachbrettfunktion» s zu 
Hilfe, welche dadurch charakterisiert ist, daB sie in allen Gitterpunkten des 
quadratischen Netzes Werte + 1 hat, und zwar so, daB die Werte in zwei be- 
nachbarten Gitterpunkten verschiedenes Vorzeichen haben. Ist nun »« Eigen- 
funktion zum Eigenwert y;, soist y= s x Eigenfunktion zum Eigenwert (8 — 4,). 
In der Tat gilt zum Beispiel im Punkt Nr. 9 von Figur 1 


4X — X49 — %q — Xg — X yg = [li Xp 
und somit 

4% + Xo + Xa + Xe t+ Hg = (8 — Ue) Xo 
oder 


4 Yo — Vio — Ya — Va — Via = (8 — Mi) Yo - 


| Speziell sind also der héchste und der niedrigste Eigenwert durch die Relation 
verkniipft 4, + “4, = 8, und der obige Fall 3 entspricht der Wahl 


je Ha Es =4, oder A= a 
Die Regel (43) fiir die Korrektur einer Versuchsfunktion v ergibt dann 
Av =— = 1, 


1) Dies pedeutet in der Figur 5, daB der Relaxationsweg in der Richtung der Ellipsenhaupt- 
achse in den Lésungspunkt einmiindet. j 

2) Vergleiche auch S. P. FRANKEL, Convergence Rates of Iterative Treatments of Partial Differen- 
tial Equations, Mathematical Tables and other Aids to Computation 4, 30 (1950). 
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also etwa im Punkt Nr. 9 der Figur 1 


1 
Avs = — r (4 Vg — Vyp — Vg — Ug Uy) - 


Der neue Funktionswert in diesem Punkt betragt also 


' 1 : 
2 Se (Vyp + Ug + Ug + Vy4) (51) 
und im Punkt Nr. 10 
aes ne 
Uio-== ey (Us + U9 + 45 + fr) » (52) 


wobei fr den gegebenen Randwert bedeutet. Dies ist das bekannte Verfahren 
von Richardson. Es ist nicht mit dem Einzelschrittverfahren von Abschnitt 1 zu 
verwechseln. Bei diesem wird die Korrektur (51) in (52) schon beriicksichtigt ; 
man hat also 5 
/ 1 Uy 

Che A (Us + Ug + 45 + FR) - 

Beim Dirichlet-Problem ist immer uw, < y,,/2, und mit wachsender Zahl der Git- 
terpunkte nahert sich uw, immer mehr dem Wert Null. Die Konvergenz des 
Richardsonschen Verfahrens ist daher hoffnungslos schlecht, da es gegen den 
SchluB hin einen ausgepragten Kafig liefert, zu dessen Beseitigung kein ent- 
scheidender Schritt getan wird. 


b) Pauschalkorrekturen 


Nach (43) und (45) erhalt man fiir den Versuchsvektor v, nach k Schritten 
durch Summation von geometrischen Reihen 


(1 —A yy)*¥—1 (1 —Au,)*—-1 


Up —- Vo ote oa xy + wae + Ky Xn e (53) 
My hen 
Die Lésung des Problems ergibt sich fiir k + co als 
x x ‘ 
SR a *s cea aca ri Xn (54) 
und damit der Fehler 
1—Aun,)" {ee ee 
Up, -—U= Oy ( r My Ny eee Oe, ( i) Xn (55) 
1 n 


Wir wollen nun annehmen, daB wir am SchluB eines Gesamtschrittverfahrens 
arbeiten und uns im Kafig befinden, so daB naherungsweise die héheren Eigen- 
funktionen in (53) und (55) weggelassen werden kénnen. Man kann nun folgen- 
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des Verfahren zur Beseitigung des Kafigs verwenden: Man bilde fiir zwei In- 
dizes 1 < k den Vektor v; — v,. Sowohl er als auch der Vektor Vv, — u sind an- 
nahernd proportional zum ersten Eigenvektor “,, also gilt angenahert 


UE 07 + C¢ (v, a Vp) 


mit einer beliebigen Konstanten c. Nimmt man daher wv; — v, als Gewichts- 
vektor # fiir eine allgemeine Relaxation im Sinne von Abschnitt 2, so wird man 
in die Nahe der Lésung kommen. Es ist ratsam, tiberhaupt jede Relaxations- 
rechnung durch eine solche Pauschalkorrektur abzuschlieBen. 


5. Das n-Schritt-Verfahren!) 


Beim Verfahren des starksten Abstiegs ist die Wahl des Gewichtsvektors 
durch den gerade vorliegenden Stand der Residuen bestimmt; sie ist also unab- 
hangig von der Vorgeschichte der ganzen Relaxationsrechnung. Beniitzt man 
diese Vorgeschichte — auch die eben hergeleitete Methode der Pauschalkorrektur 
tut dies —, so kann man ein Verfahren finden, das in endlich vielen Schritten 
zur Losung fithrt, also zugleich eine sukzessive Approximation an die Lésung 
und einen Algorithmus zu ihrer Bestimmung liefert. Es seien nun wieder vp, 
V;, ---, Vg, .-. die sukzessiven Versuchsvektoren, 7, 7,, ..., 7,, -.. die zugehérigen 
Residuenvektoren und 4p, f;, ---, Py, -.. die Gewichtsvektoren. Der Ubergang 
von VU; ZU v;,,. ist ein allgemeiner Relaxationsschritt mit dem Gewichtsvektor 
im Sinne von Abschnitt 2. Es gilt also nach (25) 


Vesa = Uz a he Pi mit An ae ae (56) 


und 
Vins = 14 Ay DP, 5 (57) 
ferner nach Satz 2 


(Feeds Px) =0. (58) 


Das Verfahren des starksten Abstiegs besteht in der Wahl £, = —7, fiir den 
jeweiligen Gewichtsvektor. Jetzt wahlen wir zwar noch 


Pe——=2 55 (59) 


aber 
Pr = — 1, + Ex Pra tip, k= 0., (60) 


1) Zusatz bei der Korrektur: Nach Niederschrift der vorliegenden Arbeit stellte ich anlaBlich 
eines Besuches im Institute for Numerical Analysis (University of California) fest, da® dieses Ver- 
fahren etwas spater auch von einer Gruppe des genannten Instituts entwickelt worden war. Ein 
vorlaufiger interner Bericht fiir das National Bureau of Standards wurde von M. R, HEstENEsS im 
August 1951 gegeben (N.A.M.L. Report, 52-9). 
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Wir fiigen also eine Korrektur hinzu, die proportional ist zum vorhergehenden 
Gewichtsvektor. ¢,_, wird so bestimmt, daB #, und p,_, konjugiert sind?) im — 
Sinne von Abschnitt 3, Untertitel b. Also 

(Px, Doxa) = 9 (61) 


oder 


( f iC? D c— ) 
—(%z, Dpy_a) + Ena (bra DPra) =9, Sea ee = (62) 


Nun folgt zunachst aus (57) unter Beriicksichtigung von (58) und (61) 


(Trias Pra) = (Yer Pea) + Ax (Doe Pia) = 9- (63) 
Indem man 7; aus (60) ausrechnet und (58), (63) verwendet: 
(Tests Ve) = —(esa Pe) + Exar (Teva, Pra) = 0- (64) 


Fir die folgenden Betrachtungen kénnen wir A; + 0 annehmen. Anderenfalls 
ist ndmlich nach (57) 7,,, = 7,, und aus (64) folgt (7;, 7) = 0, und die Relaxation 
ist beendet. Jetzt beweisen wir 

Satz 5: Beim n-Schnitt-V erfahren bilden die Gewichtsvektoren ein konjugiertes 
System. 

Beweis durch Induktion. Es sei also bereits gezeigt, daB po, py, -.., Pe 
paarweise konjugiert sind. Wir miissen zeigen, daB #,,, zu diesen Vektoren 
konjugiert ist. Wir wissen bereits ();,,, Dp,) = 0, berechnen also. zunichst 
(Per, DPp-a). Aus (57) folgt 


(Yiesas Mea) = (hes Yea) + Ax (DP x, Vea) = Ax (DP x» Te-1) (65) 
wegen (64). Fiir # = 1 ist dies 0 wegen (59); fiir k > 1 folgt aus (60) 
V1 = —Pr-1 + Exe Pr-s - 
Daher aus (65) ebenfalls 
(Teas Tra) = —Ax (DP, Pia) + An Ex-g (Dix Pr) = 0 (66) 


nach Induktionsvoraussetzung. 
Jetzt folgt aus (60) 


(Pry Dp) = — iva, DP x1) + €x (Pr Dp;,4) = — (Tava DP) (67) 


wegen der Induktionsvoraussetzung. Aus (57): 


1 
Dbpay Se — (Fe Tey) : 
sae 
1) In Figur 5 bedeutet dies, daB man vom Naherungspunkt v aus direkt auf das Zentrum der 
Ellipse hinzielt. Denn der vorhergehende Gewichtsvektor ist Tangente an die Ellipse in v. 
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Daher ergibt (67): 


(Presa Dpyz4) = — oe 


k-1 


i 
(Test Yr) ai a (Vesa Yp-1) = 0) (68) 
k-1 
wegen (64) und (66). Es sind also #,,, und #,_, konjugiert. Genau analog 
beweist man zunachst im weiteren 


(Yes %e-2) = 9 unddann (px41, Dbz) = 0 (69) 


und so fortfahrend schlieBlich den Satz 5. Aus (64), (66), (69) usw. folgt noch 

Satz 6: Bem n-Schritt-V erfahren bilden die Residuenvektoren ein orthogonales 
System. 

Ist ~ wieder die Anzahl der Unbekannten in den urspriinglich gegebenen 
_ Gleichungen (1), so kénnen in diesem Orthogonalsystem héchstens 7 von Null 
verschiedene Vektoren sein. Daher gilt 

Satz 7: Das n-Schritt-Verfahren fiihrt nach hichstens n Schritten zur exakten 
Losung der gegebenen Gleichungen. 

AuBer den in den Satzen 5, 6 genannten Orthogonalitatseigenschaften gibt 
es noch mehrere andere. Wir erwahnen nur 

Satz &: Der Residuenvektor 1;,,, ist honjugiert 2U 1, 4, ---) Vp—a 

Dies ergibt sich sofort, wenn man die 7; aus (60) berechnet und Satz 5 an- 
wendet. AuBerdem ergeben sich manche Identitaten zwischen skalaren Pro- 
dukten, zum Beispiel aus (60) 


. 


(Yn» Pu) = — (us Te) + Ena (Yes Pea) = — (ue, 1x) (70) 


wegen (58). Somit kann man an Stelle von (56) die Formel 


_ Ventre) ~ 6 fi 
Ar (DPx, Pr) 2 ( 


_benutzen. Ferner aus (57) und Satz 6: 
(Tvs Vet) = Vesa 4) + An (Ter, DPx) = An (Yer1 DP) - (72) 
Damit kann man (62) schreiben 


1 ay V1) 


eel ib Dea 
und indem man (71) einsetzt: 


Wiese) toms (73) 


Bios So 
thes (Ge Vee) 
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ZusammengefaBt besteht die Methode in folgendem Rechenzyklus: Gegeben 
seien die Residuenvektoren 7;,_,, 7, und der Gewichtsvektor p;,_,. 
1. Berechnung von &;_, aus (73) . 
2. Pe = —e + Ex-s Pes, Probe (py, Dox) = 0- 


(74) 
3. Berechnung von D#; und A, aus (70) . 


4. Vea =e + Ap Dpy, Probe (rp41, 7) = 0. 


Das Verfahren startet mit k = 0 beim Punkt 2, wobei #;_, gleich Null zu setzen 
ist. Ist dieser Zyklus oft genug durchlaufen, so summiert man nachtraglich alle 
Korrekturen Av; = Ay py. 

Wir besprechen noch zwei weitere Anwendungen des Verfahrens: 

a) Nach Satz 5 bekommen wir gratis ein System konjugierter Gewichtsvek- 
toren mitgeliefert. Wir wollen einmal annehmen, daB dieses System vollstandig 
sei, also aus 2 Vektoren bestehe. Miissen nun die Gleichungen (1) noch einmal 
aufgelést werden mit anderen absoluten Gliedern /;, so kann dies nun mit einem 
Schlage geschehen. Man setze den Lésungsvektor an 


hb = hy Po + 0 Pa + 2+ + On Pn - 
Skalare Multiplikation mit D#,; liefert dann sofort 

(u, Dp;) = (Du, pi) = % (i, Dp) , 
und aus den gegebenen Gleichungen (1’) 


(1, Pi) 
(Pe DP) oi 
Dies bedeutet auch, daB eine Simultanrelaxation im Sinne von Abschnitt 3 mit 
den Gewichten #; ausgefiihrt wurde, ausgehend von einem ersten Versuchsvek- 
tor, der identisch Null ist. Diese Simultanrelaxation kann man auch ausfiihren, 
wenn das System der Gewichtsvektoren nicht vollstandig ist. Man erhalt dann 
nur eine angenaherte Lésung. Verbessert man diese wieder nach dem 2-Schritt- 
Verfahren, so erhalt man automatisch eine Fortsetzung des Systems konju- 
gierter Gewichtsvektoren!). 
b) Es soll noch das Eigenwertproblem Dx = yu x gelést werden. Hier ver- 
wenden wir das System der Residuenvektoren, welches wir als vollstandig an- 
nehmen; wir setzen also an 


He Ue al e 
K= Colo t Sits 2° + Sala 


1) Natiirlich kann man auf diese Weise auch Matrizen invertieren. 
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De = 4% erzibt 

c v oe iS = i Ec 

Tt) Dry te ae Sn-1 Di ay = [4 (E> 0 Faves Sn Vora) : 


Durch skalare Multiplikation mit den einzelnen Residuenvektoren entsteht 
wegen Satz 8 das gestaffelte symmetrische Gleichungssystem 


(Dr. => ie) Ene =te (Dr ndis ah bnet = UL (ia ¥n—1) oper ? 


dessen Eigenwerte « berechnet werden miissen, was durch Kettenbruchent- 
wicklung in bekannter Weise leicht geschehen kann). Ahnliche Orthogonal- 
systeme, wie sie hier durch die Residuenvektoren geliefert werden, verwendet 
Lanczos?) zur Lésung von Eigenwertprobiemen. 

Die nachfolgende Tabelle enthalt ein numerisches Beispiel. In einem quadra- 
tischen Netz von 3 - 3 Gitterpunkten soll eine Gitterfunktion « so bestimmt 
werden, daB 


Du = 10000 in der linken oberen Ecke, 
Dit =O in den iibrigen 8 Gitterpunkten. 


Dabei ist D der Dirichlet-Operator von Figur 2. Es handelt sich also um die 
Bestimmung einer Greenschen Funktion. Es wurde mit der Versuchsfunktion 
Vp = 0 gestartet. Jedes Kastchen der Tabelle entspricht einem Gitterpunkt, 
und es sind in ihm die jeweiligen Werte von f,, Dp, 7,41 in dieser Reihenfolge 
auf einer Zeile protokolliert. Die Rechnung schreitet nach unten fort. Am 
SchluB sind die Werte der gesuchten Funktion « angegeben. Infolge der spe- 
ziellen Symmetrieverhaltnisse braucht das Verfahren nur sechs Schritte statt 
deren neun. AuBer den in (74) erwahnten quadratischen Proben kann man noch 
lineare Proben vom Typus (15) benutzen. 

Das n-Schritt-Verfahren ist besonders bequem, wenn es sich um einen Dz/- 
ferenzenoperator D handelt, weil dann die Berechnung von Dp viel weniger Mithe 
macht, als wenn ein ausgefiilltes Gleichungssystem (1) vorgelegt ist. 


1) L. Cotzatz, Eigenwertprobleme und ihre numerische Behandlung (Akademische Verlagsgesell- 


schaft, Leipzig 1945), S. 316. 
2) C. Lanczos, An Iteration Method for the Solution of the Eigenvalue Problem of Linear Differ- 


ential and Integral Operators, J. Res. Nat. Bureau of Standards 45, Nr. 4 (1950). 


28 EDUARD STIEFEL ZAMP 


k | Pre | Dex Th+1 Pr: | mn | e+ Pr Dey | "e+a 
— 10000 0 

0 | +10000 | +40000 0 0} —10000 | — 2500 0 ) 
1 | + 1250 0 0 | +2500] + 8750 ) 0 | —2500] —714 
2.\ + 306 0 O}+ 612] — 1 o| +714 
ek} 0 0}+ 186 0 0 | +217 
4|/+ 20 ) Ol+ 40} = 1 O| + 47 

| 2 | | | 982 | | 312 

0 0 

0 0 | —10000 | — 2500 ) 0 0 
1 | + 2500 | + 8750 0 0} — 5000 | —1429 0 0 ) 
2\4 612 | — 1 0 | +1429] + 4492 0 0.| —2143 | —682 
3 |e 186 0 O}+ 434 0 0 | +682 | +2077 0 
Gal ee Oe ee 1 Os] Ot 1 2 8 AS CY Ie pee ) 

| 982 | 
0 0 
i Gee 
Bie 9714 4.2 
3|/+ 217 —1364| — 448 
4 0| +448 | +1498 0 


or 


(Po, Dp) = 400-108 A, =0,25, &) = 0,125, 
(ry, %4) =12,5-10°, (p,, Dp.) = 4375-108 A,=0,28571, 2, = 0,244.94, 
(7s, ¥0) = 3061633, (fp, Dp») = 9622276, 4, — 0,31818, ¢, = 0,30384, 
(rs, 7s) = 930248,  (p,, Dps) = 2833028, A, = 0,32836, 2, — 0,21575. 
(74,74) = 200704, (py, Dp) = 671012, A,—0,29911, 


Yo, %) = 100 - 108, 


Ein groBeres Beispiel dieser Art, das auf 106 lineare Gleichungen fiihrte, 
wurde im Institut fiir angewandte Mathematik der Eidgendéssischen Techni- 
schen Hochschule mit Hilfe der programmgesteuerten Zuse-Rechenmaschine 
durchgeftihrt. Es handelt sich um folgendes physikalisches Problem: 
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Auf eine parallelogrammférmige Scheibe, die durch ein Gelenk A und eine 
Pendelstiitze B gehalten wird, wirkt eine Einzelkraft P in ihrer Mittelebene. 
_ Man sucht die Spannungen und Verschiebungen. 

Um wenigstens eine der numerisch ungiinstigen Differentiationen zu ver- 
meiden, berechneten wir statt der Airyschen Spannungsfunktion direkt die 
_ Verschiebungen, indem die diesbeziigliche Randwertaufgabe fiir ein System 


von zwei partiellen Differentialgleichungen mit Hilfe der Differenzenrechnung 
approximativ gelést wurde. Die Wahl des eingezeichneten Dreiecksnetzes mit 
53 Gitterpunkten (der festgehaltene Punkt A kann weggelassen werden) fiihrt 
so auf 106 lineare Gleichungen vom Typus (1) fiir die Verschiebungskompo- 
nenten in Richtung der Koordinatenachsen in den Gitterpunkten. Fir den 
Start des u-Schritt-Verfahrens wurden die Verschiebungen als Null angenom- 
men und anschlieBend 90 Rechenzyklen durchgefithrt. Nach dem neunzigsten 
Schritt waren die Residuen so klein, daB die Rundungsfehler zu iberwiegen 
begannen und somit abgebrochen werden konnte. (Die SchluBresiduen waren 
ungefahr der 10-te Teil der Ausgangsresiduen.) Infolge der beschrankten Spei- 
cherkapazitat der Maschine dauerte ein Zyklus etwa 2h 20m. Die Durch- 


} 
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rechnung eines zweiten Lastfalles gema8 Formel (74) ergab nur eine angenad-— 
herte Lésung, die durch weitere Relaxation verbessert werden muBte. Dies ist 
darauf zuriickzufiihren, daB das konjugierte Funktionensystem der #; infolge 
der nur fiinfstelligen Rechnung zu stark mit Rundungsfehlern behaftet war. 


6. Harmonische Gewichtsfunktionen 


Es soll hier am speziellen Beispiel des Dirichlet-Problems eine Querverbin- 
dung zu den Methoden von TREFFTZ, BERGMANN und Picone?) hergestellt 
werden, welche darin bestehen, Partikularlésungen der gegebenen partiellen 


Differentialgleichungen zur Lésung von Randwertaufgaben heranzuziehen. Wir — 
betrachten zu diesem Zweck in einem unendlichen, die ganze Ebene ausfiillenden _ 


quadratischen Netz Gitterfunktionen #, die in jedem Punkt die Differenzen- 
gleichung Dp = 0 erfiillen. Dabei ist D der Operator von Figur 2a. Wir nennen 
dann # eine «differenzharmonische» Funktion. So sind zum Beispiel Real- und 
Imaginarteil des Ausdrucks 


m+n=k 


fiir k = 0,1, 2,... solche differenzharmonische Gewichtsfunktionen2), die den 
harmonischen Polynomen bei der Laplaceschen Differentialgleichung entspre- 
chen. (Netz und Koordinatensystem wie in Figur 8; m, ganze Zahlen = 0, 
? = imaginare Einheit). 

Ist nun etwa unser altes Dirichlet-Problem von Figur 1 vorgelegt, so kann 
man, ausgehend von der Versuchsfunktion v= 0 als Gewichtsfunktion fiir die 
Relaxation, die Werte einer differenzharmonischen Funktion in den 25 inneren 
Gitterpunkten nehmen. (Am Rande des Grundgebiets und auBerhalb desselben 
ist die Gewichtsfunktion natiirlich gleich Null zu setzen.) 

Relaxiert man ausschlieBlich mit differenzharmonischen Gewichtsfunktio- 
nen, so kénnen nur Residuen in den Punkten vom Typus Nr. 10 und Nr. 5, also 
entlang des Randes auftreten; der innere Kern von 3 - 3 Gitterpunkten bleibt 
dauernd frei von Residuen. Dementsprechend braucht man im Beispiel der 
Figur 1 genau 16 linear unabhangige differenzharmonische Gewichte, um das 
Problem zu lésen. Mit anderen Worten: Im urspriinglich gegebenen Problem 
ist die Anzahl der Unbekannten gleich der Zahl der inneren Punkte im Grund- 


1) E. Trerrrz, K onvergens und Fehlerabschatzung beim Ritzschen V erfahren, Math. Ann. 1 00, 503. 
(1928). — S. BERGMANN, Punch-card Methods applied to the Solution of the Torsion Problem, Quart. 
appl. Math. 5, 69 (1947). — S. BERGMANN und M. Scuirrer, Kernel Functions in the Theory of 
Partial Differential Equations of the Elliptic Type, Duke Math. J. 15, 585 (1948). — M. PICONE, 
Appunti di Analisi superiore, Rondinelle 1940. — G. FIcHERA, On some general Integration Methods, 
ted tm Connection with Linear Differential Equations, Pubbl. Ist. naz. Appl. Calcolo, Nr. 245. 
1950). 

*) Mitteilung von Herrn Dr. H. RUTISHAUSER. 


ee 
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gebiet (hier 25) ; die eben entwickelte Relaxationsmethode bewirkt eine Reduk- 
tion auf 16 Unbekannte (= Zahl der «randnahen» Punkte). 

Man kann mit differenzharmonischen Funktionen als Gewichten einzeln und 
simultan relaxieren und schlieBlich auch durch den Schwarzschen Orthogonali- 
sierungsprozeB ein System von konjugierten differenzharmonischen Gewichts- 
funktionen herstellen. Diese Methode entspricht der Bildung der Bergmannschen 
Kernfunktion fiir die Laplacesche Differentialgleichung. 

Die Ausgestaltung dieses Ansatzes, der auch im Zusammenhang mit der 
nach FRIEDRICHS benannten Transformation) von Variationsproblemen inter- 
essant ist, muB weiteren Arbeiten vorbehalten bleiben. Dasselbe gilt von der 
Verwendung von Grundlisungen als Gewichtsfunktionen. (Unter einer Grund- 
losung verstehen wir eine im unendlichen Gitter definierte Funktion #, die in 
allen Gitterpunkten mit Ausnahme eines einzigen «Quellpunktes» differenz- 
harmonisch ist. Im Quellpunkt soll Dp = 1 sein, wobei D wieder der Operator 
von Figur 2a ist.) 


7. Fehlerabschatzungen 


Es soll aus der Gr6éBe der Residuen in einer Relaxationsrechnung auf die 
GroBe der Fehler in den Unbekannten geschlossen werden, um festzustellen, 
wann die Relaxation abgebrochen werden kann. Es sei wieder wu der Lésungs- 
vektor, v der Versuchsvektor und 7 der zugehorige Residuenvektor, so daB nach 
(1’), (3') gilt 

Du+l=0, Dv+l=r. (76) 


Fiir den Fehlervektor f = v — u hat man also 
Di=r. (77) 
Wir fiihren nun die Schwarzschen Konstanten?) von / ein: 


=), 
=(/, D) =(6.7), 
= (f, D*f) = (DL DN = (7), 7 Bs 
= (f, D®#) = (Df, D®f)= (7, Dn), 
n=, Df) = (D*f, Df) = (Dr, Dn), 


el bids) (e: inei shes sive) (4.6) (ele! ©! alle. 'm egieniel «like -e] (ela) eelione re!) e) (eis 


1) Vel. R. Courant und D. Hirsertr, Methoden der mathematischen Physik, Band 1, 2. Aufl. 


Springer, Berlin 1931). 
& aol L. Coriarz, Eigenwertprobleme und thre numerische Behandlung (Akademische Ver- 


lagsgesellschaft Leipzig, 1945), 4. Kapitel, § 12. 
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Ferner die Schwarzschen Quotienten 


ee Diet (k= 0, 1,2, 028) (79) 


1: 


welche bekanntlich eine monoton wachsende Folge bilden. Das geometrische 


Mittel zwischen o, und o, betragt 
sao Von See. “ 
Fiihren wir noch den mittleren quadratischen Fehler 
ye (eee 


und das mittlere quadratische Residuum 9 = V(r, In ein, so gilt also 


oO 


bis ate (81) 


O7i/2 


Auf dieser Formel beruhen alle Fehlerabschatzungen. So folgt zum Beispiel 
aus der erwahnten Monotonitat 


O12 SO, So, = (in Riv) (82) 


Q 
Go RO) * (83) 
Ferner 


wenn /1, den kleinsten Eigenwert des Problems Dx = wx oder der quadratischen 
Form Q(x) bedeutet. Somit 


i (84) 


Aus (83) lernen wir folgende wichtige 

Regel: Um eine Relaxationsrechnung giinstig abzuschlieBen, muf die ver- 
bleibende Residuenverteilung einen miglichst hohen Rayleighschen Quotienten 
haben, das hei ft wie der hichste Eigenvektor aussehen. 

Beim Dirichlet-Problem bedeutet das, daB die Vorzeichen der verbleibenden 
Residuen schachbrettartig abwechseln sollten (vgl. Abschnitt 4). Hért man mit 
einem Kafig auf, iiberwiegt also etwa das positive Vorzeichen in der Residuen- 
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verteilung, so kénnen grobe Irrtiimer entstehen, indem dann nach (83) 6 um 
einige Zehnerpotenzen gréBer sein kann als g. Ahnliches gilt fiir die biharmoni- 
schen Probleme der Airyschen Spannungsfunktion. 

Die Abschatzung (84) ist in der Praxis meistens unbrauchbar. Erstens weiB 
man bei komplizierten Problemen zu wenig tiber j,, und zweitens beriicksich- 
tigt (84) nicht die Gestalt der Residuenverteilung, ist also auf den schlimmsten 
Fall des ausgepragten Kafigs zugeschnitten und gibt daher meistens ein ent- 
mutigendes Resultat. 

Um eine etwas bessere Information als (83) iiber 6 zu erhalten, kann man 
folgendermaBen vorgehen. Man bildet auBer o, = R(v) auch noch 


(Dr, Dr) 


(7, Dr) 


03 = 


und berechnet dann angenahert oj). durch logarithmische Extrapolation aus 
6, und o3: 


| w 


log Oy. = logo, — = (logo, — loga,) . (85) 


Auch mit dieser Formel wird der Fehler meistens noch unterschatzt. Hingegen 
gilt in (83) das Gleichheitszeichen, wenn 7 irgendein Eigenvektor ist, weil dann 
die Schwarzschen Quotienten gleich dem zugehoérigen Eigenwert sind. Beim 
n-Schritt-Verfahren von Abschnitt 5 kann man die Schwarzschen Quotienten 
allein aus den Zahlen A,, €, berechnen. 


Summary 


The general outlines of the so-called relaxation technique are developed. By 
“relaxation” is meant every ‘‘step-by-step procedure’’ for solving systems of 
linear equations based on minimizing quadratic forms. After a short discussion of 
the trial methods developed by Southwell and his school, allowing full leeway to 
the intuition of the computing person, the general mathematical background is 
treated. § 3, 4 are the central parts of the paper. After a study of the gradient 
method it is shown that relaxation methods are not necessarily successive approx- 
imations taking an infinite number of steps but that it is possible to speed up 
convergence such that the desired result is reached in a finite number of steps. 
These methods may be suitable for use on sequence-controlled computing machines. 
Special consideration is given to the well-known fact that relaxation very quickly 
smoothes the trial function but that it may be a combersome task to get rid of 
the remaining smooth residual distribution. 


(Eingegangen: 14. 7. 1951.) 
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4 ee 
Unsteady Flows with Free Boundaries’) 
By Davip GiLBarG, Bloomington, Ind., U.S.A.*) 


Unsteady flows with free boundaries are of interest in the theory of under- 
water cavity formation. Mathematically they present a generalization of the 
now classical Helmholtz problem of conformal mapping. The time dependent 
character of the flows introduces two difficulties beyond those already apparent 
in the stationary flow problem; namely, (1) the unsteady form of BERNOULLI’s 
equation replaces the much simpler condition of constant flow speed as the 
kinematic boundary condition on the constant pressure surface (free bound-_ 
ary); and (2) the free boundary, being a material line, is generally not a stream- 
line. The effect of these difficulties is to place the exact solution of the unsteady 
free boundary problem beyond the realm of present techniques. However, as 
is shown here, we can achieve considerable mathematical simplification, with- 
out doing great violence to the physical problem, by replacing the requirement 
that the free boundary be a material line with the approximating condition 
that it be a streamline. It is physically reasonable, at least for not too rapidly 
varying flows, that this model of the unsteady cavity should yield results in 
substantial agreement with the exact theory. We shall see that the unsteady 
flow problem, thus modified, is not essentially more difficult than the station- 
ary free boundary problem’), and, as regards flows past polygonal obstacles, 
that it can be solved completely by standard conformal mapping techniques. 

Von KArmAN‘) has exhibited an explicit example of a symmetric unsteady 
flow with finite closed wake (or cavity) behind a flat plate. We intend here to 
solve in detail (subject to the above modification) the problem proposed by 
him of determining all symmetric unsteady flows with finite closed cavity 
behind a flat plate, and more generally, behind any symmetrical polygon. 
Among the solutions thus obtained will be a class of flows with constant cavity 
shape (of which von KARMAN’s example is a particular case), for which the 
free boundaries are indeed simultaneously material lines and streamlines, and 
for which, therefore, the present model and the exact theory coincide. The 
solution of the asymmetric flow problem, which is more difficult only in detail 
and not in principle, is briefly indicated here. 


1) Prepared under Navy Contract with Indiana University. 
*) Graduate Institute for Applied Mathematics, Indiana University. 


3) For this see, e.g., L. M. Mirne-Tuomson, Theoretical Hydrodynamics (Macmillan, New 
York 1938), Chapter 12. 


4) Von KArmAn, Ann. Mat. pura appl. [4] 29, 247-249 (1949). 
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Finite cavity flow past a flat plate 


Consider the problem of symmetric finite closed cavity behind a stationary 
flat plate. Let the flow be described by a complex velocity potential, 


w(z, t) = p(x, y,t) +7 p(x, y, 8), (2=%-+7 y, ¢= time), 


which is assumed to have the properties: 

(Gi) lim wv’ = U(d), (w’ = 0w/dz); 

(ii) w’ + 0 in the interior of the flow region (i.e., there are no stagnation 
points in the interior of the flow) ; 

(ui) free boundaries, on which pressure is constant, separate from the ends 
of the flat plate and meet to enclose with the plate a finite region (the cavity) 
(figure 1a); the potential w is defined in the entire z plane outside this region; 

(iv) the free boundaries are streamlines, namely curves on which w(x, y, #) 
= const. 

On the free streamlines, BERNOULLI’S equation takes the form 


SP 4 > \w! | = olf). (1) 


In most applications, c(t) = U?/2+ p..— p,, where p, and #, are the static 
and cavity pressures, respectively (assumed to be either constant or known 
functions of time only). Equation (1) for nonstationary flows replaces the 
condition, | w’|=const on free streamlines, for the classical steady flow problem. 

Because of the postulated symmetry we need consider only the upper half 
of the flow plane, and corresponding to it, the upper half-plane of the complex 


Fig. la 
Finite cavity behind a flat plate. 


potential, w. If J represents the point at infinity in the physical plane, A the 
stagnation point on the plate, B the upper edge of the plate, and C the terminus 
of the free streamlines (figure 1a), then the dividing streamline 1A BCI maps 
into the real (~)-axis of the w-plane (figure 1b); the points wy, Wg, We, vary 


: 
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zero for all time. 


with ¢ in general, although any one of them, say wz, may be taken equal to 
Let now : 


weal [1->(¢+=)]- (> 0) (2), 


define a mapping of the upper half unit circle, |¢| <1, Im¢ = 0 (figure Ley 


: 
, 
t-plane — 
: 
: 
; 
| 
: 
-| +1 7 
(C) (1) (A) (B) 
Fig. 10 Fig. le ; : 
Plane of the complex velocity potential w. Plane of the auxiliary variable C. : 


into the closed upper half w-plane, Imw = 0. Under this mapping we have 
the correspondence, for a suitable value of « at each ¢, 


= 


$=0++ 0,00, wp=0->f=1, wef =-1, 
w,<> €=hk (= areal number, 0 < k< 1); 


the free streamline BC is the image under (2) of the semicircle, |¢| = 1. 
We define, | 
w(C) = logw’ = log|w’| — 7 6, 
where @ is the inclination of the streamline through the point z at time ¢. On 


the real ¢-axis, 9=0 for -1S0< k, 0=2/2 for k<¢ <1. Hence the 
function, 


2 x k— ¢ \1/2 
Ie) = o(f) — log (°°), 


is analytic in the ¢ half-circle (recall w’ + 0), and real and continuous on the 
real axis. By Schwarzian reflection /(¢) can therefore be defined in the full 
unit circle |¢| <1. On |¢| = 1, we have by (1), 


Re /(C) = Re w(C) = log |w’| 


=> log{—2 a'() E ~ Stes =)| + 2 o)} (3) 
= + log] 4@ + BY) (¢+=)}, (]¢|=1) 
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ai 


where 
A(t) = —2a'(t)+2c(t), Bt) =a'(t). 


Since f(¢) for |¢| <1 is uniquely determined by its real part on |¢| = 1, it 
suffices to find an analytic function g(f) with the properties: (a) Le() |=. 
on |¢| = 1; (b) h(C) = g(¢) [AC + B (2 + 1)] is regular for |¢| <1, without 
zeros or poles; (c) A(¢) is positive for ¢ real. For, given such a function g(C), 
it follows by (3) that log h(¢)/2 has the same real part on |¢| = 1 as f(), is 
regular for |€| < 1 and real for ¢ real, and hence 


Ho) = > log felt) [46+ B(C?+ 1)]}. 


The function g(¢) is easily determined. Let €) be a zero of the quadratic 
B(¢?+ 1) + A(¢) = 0, with |¢,| <1 (the case B = 0 is obviously degenerate 
and of no interest). We consider the two possibilities, 


Ticiel Sool OL iain 


Case I: If |€)| < 1, then ¢) = g, a real number; the function g(¢) satisfying 
conditions (a) and (b) above must map the unit circle into its exterior and 
have a simple pole at ¢ = ¢, = 9; thus g(¢) is of the form e4(o € — 1)/(f — @). 
To satisfy (c), h(¢) = B(E —1/e) (oe € — 1) e” is real and positive for ¢ real if 
and only if Be'*/o is real and positive; hence e** = sign (B/o), and thus 


fle) = = log[|F| A — 0"): (4) 


seal lait) Calder Cone’? Mandisincesa(C) BC <5) (C — C,) e(G), it 
follows that g(£) = sign B satisfies conditions (a) to (c), and 


HC) = + log {| B| (C2 2 cos B + 1}. (5) 


From the above, we have, therefore, in case I, 


, | Bir eae 


and, since w’(0) = U, then| B| = U?|o|/k, and 


1— ¢/k \1/2 
(Case There Ul oe) Ge a (6) 
Similarly, in case II,|B| = U*/k, and 


Se Oe 
(Case 11) w'= U(f?—2Ccosp+yP(F——ey. 0) 
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To determine k and «, and from these the entire flow, it remains only to’ 
satisfy the conditions that z¢ lie on the real axis, and to fix the mappings so 
that zp = hi, h being the half-width of the plate. The former condition requires 
that z(C) be single-valued in the neighborhood of ¢ = 0 (z = co in the physical 
plane), or, in other words, that the residue of dz/df = 0 at € = 0. From (2), 
(6), (7) one obtains readily, 

(Case I) k ~~ = Zo, (Case ll)"k aes = "2 CSD: 
In both cases, we have therefore —1< g, cos B < 0, and ¥2—1 hoa 
The remaining condition states, 


1 


; US UD, 


k 
oe ] op ee a 9 ge 
(Case I) = Co) f = a3 5 ( aE =" df= ta) F(R). 
: e-F0-F] 
(et) = (28) frm en nee ee 


1 1 . 

[cos B = > (2 -;)| (9) 
Both 44(k) and 4,(k) are expressible in terms of elementary functions. The . 
unknown parameter k is now determined within a constant (its initial value) — 
in terms of U(f) after elimination of « from the relations, . 


al hy deh. Up 3 nd 
(Cael) a= aay \Gl=|Bl=-42=5 02 (4, —1). (10) 
im _ Uh | da ys Ge 
(Case Il) a= pe, ra = (Bl. (11) 


With the determination of k and hence «, the flow is given to within a 
quadrature by (2) in conjunction with (6) or (7). 

To summarize; the symmetric unsteady flows with finite cavity behind a flat 
plate of width 2h are given in parametric form by 


Uh iis 
(Case I) : < F(R) E ~F (<¢ - =)| ; 
Ww 


* Sp La CRA sky, it 
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Uh 1 1 
| sare ed Gasca le 
(Case) ) : : ae eee eee ne 
ae (¢2— 2 cosp +1) Geaaa |cos B = + (k —=)|; 


| a Uh ye 2 il 
Bore) (i) 2 U* (a 1), 
d Uh Ui? 
(Case IT) lS FE) ; 


The equation of the free boundary is given by: 


F, 2 (e*®) 


Dee es 


(OCS S22) 


Cavities with constant shape. The unsteady flows whose free boundaries are 
of constant shape are given exactly by this model, since for such flows the free 
streamlines are also material lines. We have from the preceding that these flows 
are characterized by the condition k = const. From (10) it follows 


al 


ser [®=7R@ (Ge-1)) a 


U = —— 
Us Sei 


the sign in the denominator being opposite that of da«/dt (= B) and of dU/dt. 
The particular flow for which 9 = —1, k= V2 —1, and the sign in (12) is 
negative, is the example given by voN KArMAN!). We note, contrary to his 
conjecture, that this flow is only one of a class of flows with constant cavity. 
Similarly, such a class can be constructed based on (11). 

Shape of the cavity. From the fact that 6 = —argw’=0 at ¢=—1 if 
o + —1 in (6) we conclude that the cavity shapes defined by the above flows 
have cusps at the rear. In the instance 9 = —1 in (6) and 8 =z in (7) the 
flow (which is that given by von KArMAN), has a stagnation point in the rear 
and is of constant shape. 

For general flows of type II an examination of formula (7) shows that 
§ = —Imlogw’ is a decreasing function of s (= arg ¢) on |¢| = 1, except at 
s = B where 6 makes a jump increase of z/2; since @ has initial value z/2 at 
s = 0 and final value 0 at s =z it follows that the free streamlines must at 
some point cross the axis of symmetry of the flow. Thus the complete flow plane 
must contain a doubly covered subregion (figure 2), from which it appears 
that flows of type II are physically unrealistic. 


1) Yon KArmAn, Ann. Mat. pura appl. [4] 29, 247-249 (1949). 
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Fig. 2 


Flow of type II, showing doubly covered flow region. 
Similarly, one observes from (6) that 


nee : /k—€ \"e2 
0(s) = — Im log w leuet® == arctan 1 — 9 coss — ale log Gen iB -} , 


Simple computation shows that 6”(s) > 0 for 0 << s <a, and since @ has a 
minimum value (< 0) in 0 <s <7z, the free streamline must be an arc with 
one point of inflection, lying entirely above the axis of symmetry of the flow 
(figure 1a). This is the characteristic shape of the cavity for the flows of type I 
behind a flat plate. 


Arbitrary polygonal obstacles 


The preceding results can be extended easily to include arbitrary symmetric 
polygonal obstacles, and with some additional formal complication, also asym- 
metric polygonal obstacles. We outline these extensions briefly. 

Consider the upper half of a polygonal obstacle which is symmetric with 
respect to the x-axis, the direction of the flow. Let the angle between the axis 
of the flow and the face of the obstacle along which it divides be «%) a, and the 
vertex angles of the polygon, a; (¢ = 1, ..., 7 —1; = number of faces), 
where the angles taken are those facing the interior of the flow. Then the 
same considerations as for the flat plate lead to the following analytic descrip- 
tion of the flow: 


= eS arg 
es wei é m—-1, — tlk, 5) || eas dt Im € = 0; (13) 
= U0-e dT Goa) 
(only flows of type I are considered, for sake of brevity). The parameters 
x, k;, @ are determined as before from the conditions that the residue of 
dz/df = 0 at ¢ = 0, and that the vertices of the polygon are at assigned points 
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z; (( = 0, 1, ..., m) in the z-plane. The former condition gives readily 
mill 1 
o= S (1—a) (ke — 5). (14) 
t=1 Z 


If we introduce the function F(¢), defined by the indefinite integral 


G 
7 el ee PSO nates: 
es | CU(1— ef) Adi ae tk, ) aC, 


then the second of the above conditions gives, 
pele LR) = (Rol = 051, 25s ky = 1) > (15) 


The conditions (14) and (15), combined with the relation, 


Laer 3" De 
dt n—1 j20-a) 2 


1=0 


which is obtained in the process of derivation of (13), suffice to determine the 
quantities, «, 0, k; («= 0, 1, ..., m—1), as functions of ¢ once their initial 
values are given. Their numerical determination is of course a formidable 
problem of calculation in general. 

The problem of cavity shape is somewhat more complicated in general than 
for the simple case of the flat plate, and we forego the discussion here. 

Concerning the asymmetric flow problem, we make only the following brief 
remarks. It is now necessary to consider the entire flow plane, and its image 
in the w-plane. The latter is simply the full w-plane slit along a segment of the 
real axis; the two edges of the slit correspond to the boundary of the flow 
region, that is, obstacle plus free streamlines, each edge corresponding to one 
of the two streamline arcs into which the dividing streamline separates on the 
boundary. A mapping of the form 


a 

= TH BF feos sy — (C+ 1/0/21 (te) 
where «, f are real (functions of ¢) and 0 < s9(¢) <a, takes the upper half 
unit circle |¢| <1, Im¢é = 0, into the w-plane, the semicircle |¢| = 1 cor- 
responding to the free boundary under the mapping, and the real axis to the 
fixed boundary; the point at infinity is an interior point of the mapping. The 
mapping (16) replaces (2) in the analysis of the asymmetric flow problem. Now 
following a line of argument similar to that in the symmetric case, and deviating 
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only in the details rather than in basic method, we arrive at a (more complicated). 


formula for w’ analogous to (13), with appropriate side conditions to determine 
the parameters of the problem. However, we omit the formal details here. 

In conclusion, the author wishes to acknowledge the able assistance and 
contributions of JAMES SERRIN in this work. 


Zusammenfassung 


Instationare Stromungen mit endlich begrenztem Hohlraum hinter vieleckigen 
Hindernissen werden unter der Voraussetzung, da8& freie Rander Stromlinien sind, 
betrachtet. Das mathematische Problem wird in expliziter Form nach dem tib- 
lichen Verfahren der konformen Abbildung gelést und zeigt nicht wesentlich mehr 
Schwierigkeiten als das klassische Problem der freien Rander fiir stationare Str6é- 
mungen. Die Form des Hohlraums fiir symmetrische Str6mungen hinter einer 


ebenen Platte wird diskutiert, und Str6mungen mit konstant bleibendem Hohl-— 


raum werden mit expliziten Formeln naher umschrieben. 


(Received: April 9, 1951). 


Sur la thermodynamique des processus irréversibles 


Par KyRILLEe Poporr, Sofia?) 


Introduction 


ONSAGER®) a esquissé les fondements d’une théorie des processus irréver- 


sibles basée sur la théorie des probabilités. CAsrmir*), DE GRoor?) et d’autres | 
ont repris et développé cette théorie. Les phénoménes irréversibles s’expriment, © 


en premiére approximation, par des relations phénoménologiques du type 


J:= D7 Lig Xx; 


k 


on établit, en se basant sur la théorie des probabilités, qu’on a L,, = L,,. Pour 
arriver a ce résultat, on considére la marche des processus thermodynamiques 
dans un systéme adiabatiquement isolé, au voisinage de l’état stable. Soient 
S Ey, ty &, les paramétres définissant l’état du Systeme; 0 (69, 25-4 eee 
l’entropie du systéme 4 un moment donné ¢ et &, é%, ..., €@ les valeurs de ces 
parameétres au moment ¢= 0, correspondant a l’état stable, défini par le 


1) Institut de Mathématiques a l'Université de Sofia. 
2) L. OnsAGER, Reciprocal Relations in Irreversibl 
f ONS In stble Processes, Phys. Rev. 37, 405- g 
38, 2265-2279 (1981). re - ite 
3) H. B.C. Casimir, On Onsager’s Principle of Microscopi bili. 
, ‘ ‘ ; icroscopic Reversibility, Rev. mod ‘ 
343-350 (1945); Philips Res. Rep. 1, 185 (1946). P YW modern Phys. 17, 


4 a 
) S.R. pe Groot, Phermodynamics of Irreversible Processes (North-Holland Publishing Com- 
pany, Amsterdam, 1951). 


= 
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maximum de S. On a pour la variation AS de l’entropie 


ANS ESS ee eee see eS HAE int 
1 02S 1 (1) 
= a8, 06, GB E—B=- FS bam, | 
ou l’on a mis 
CR eS 


Les coefficients gj, = 0?S9/(0&; 0,,) désignant les dérivées secondes de S pour 


tos, ---, 62, ON a par conséquent g;, = gus 


Ti Stites 
De l’équation (1) on tire d’abord 


E 
6(AS) = = Pin hy OX 
et ensuite, en introduisant les dérivées généralisées) par rapport a ¢ 
AS = = J bux ee. (3) 
(3) peuvent étre mises sous la forme 


= DEES: OX; , 
k 


(2) 


Les équations (2) et 


avec 
yee ses) 


ONSAGER et DE GROOT appellent les dérivées généralisées x; des «flux » 


quils désignent par J;,,, les X;, — des «forces » et posent 
l= D>) Le Ay: (6) 
i 


D’aprés l’équation (4), 6(4S) peut étre interprété comme le travail des 


’ x od . 
forces X; correspondant aux déplacements respectifs 6x;,. Dans ce qui suit 
nous montrons qu’en considérant les X,, comme des forces Newtoniennes, c’est- 
a-dire en partant des équations 
2% 
aK pat (7) 
dt? 

1) Kyritte Poporr, Sur une extension de la notion de dérivée, C. rx. Acad. Sci. Paris 206, 1440 a 
1442 (1938); 207, 110-112 (1938) ; 209, 472-474, 668-670 (1939) ; Uber die verallgemeinerten A blettun- 
gen, Abh. Preu8. Akad. Wiss., math.-naturw. Klasse Nr. 2 (1942); Sur une extension de la notion de 


k»> 


dérivée, Mh. Math. Phys. 51, 115-152 (1944). 
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on arrive aux relations phénoménologiques 
Ji = ys Lix Xx 
ke 
avec 
Liz ae Ty; 


Pour arriver ace résultat, il n’est pas nécessaire de limiter ses considérations, 
au voisinage de l’état stable du systéme. Nous admettons d’une maniére plus 
générale que le moment ¢ = 0 correspond 4 un moment quelconque et que, 
par conséquent, on a, avec une approximation suffisante, adoptée par ONSAGER, 
DE GROOT et d’autres, 


1 . 
AS =— J! 8, %.— > D) ie %i Xt (8) 
Rk i,k 
avec 

Sik = Ski - (9) 
Cette généralisation s’impose par le fait que les équations d’ONSAGER et de 
DE GROOT, €étant linéaires, correspondent a la valeur infinie du temps ¢ de l’état 

d’équilibre qu’elles définissent. 

Etablissement des relations phénoménologiques 


En partant de (8), posons a l’instar d’ONSAGER et de DE GROOT 


0(AS) 
Brum re 


ce qui donne aussi dans le cas général 


(AS) =) X, Ome, ee eae 
k k 


X;, = 


Des équations 
d 2%, 


aa = Xt (7b) 


l’on tire, en multipliant par x, = «, = dx,/dt, 


et par conséquent 

oe 12 
>—— =AS+C, (11) 
C désignant une constante d’intégration. L’équation (11) exprime ainsi une loi 


analogue a la lot de conservation de V’énergie. 


Apres cette remarque préliminaire admettons, pour ne pas compliquer 
inutilement l’exposition, y= 2 et posons x,= x, %,»= y, X,=X, Ag = 
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On aura a considérer ainsi le systéme d’équations différentielles, avec une 
€criture un peu changée, 


Gx dy 


aK OY = X , Fo a BP eee (12) 


J, m, a, b, c étant des constantes, correspondant aux valeurs des dérivées pre- 
mieéres et secondes de S au moment ¢ = 0. 

De ces équations, on tire d’abord, en omettant les constantes d’intégration, 
ce qui revient a ce que les forces nulles correspondent aux flux nuls, 


Gea titaf xdted fy dt, 
t 
et, en intégrant par partie, 


ax "dy 
= (ltaxtby)t afta dt bf tap dt. 


Une nouvelle intégration par partie donne 


A 49\ ax b yg ay 2 ax A 
pres ts oe kets 5 ft Peg hoes 


En substituant sous les signes d’intégrale aux dérivées secondes de x et y 
leurs valeurs respectives X et Y et en intégrant par partie, on obtient 


dx. by, dy Ge a ak ae ees 
(1+ 4 ) te iteg me P\XES 8 Y ee ad 
sf Pg & 


En substituant aux dérivées de X et de Y leurs valeurs tirées de (12), 
aprés une nouvelle intégration par partie et des simplifications évidentes, on a 


d f 2 
(1 a8 of #2 4 a® a Os i) ae ut i b a, ou ee is) a os (¢ ae +; i) xX 


b a a? x d*y 
Beit vehias if #(0 ae te a). 


En substituant de nouveau aux dérivées secondes de x et de y leurs valeurs 
tirées de (12), on obtient 


a a2 + b? b (a a oe) dy 
(1 a eee a “) a he (ai = iran, i) BPE 
2 2 b b 
a (: hice Bate - 2 aL eee £ ? Xe de {8 ees ; 4 We (13a) 


wah fs a8 ie et Bere fs at 


et ainsi de suite, l’itération pouvant étre poussée a l’infinie. 
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Par un procédé analogue, en partant de la seconde des équations (12), on 
obtient de méme 


ate) dx ee ely: | 
" Df is et an t 


Okara b (a + c) ax b2 + ¢? dy 
(sr # + ot i) SF + (1+ ope + a i*) rf | 
b ¥ b? 
Ss (sr foes: ae a) i) Bes (¢ ry < i 4 = is u Y (13b) 


L’étude de la convergence des séries 


Nous allons montrer que les séries qui figurent dans ces développements 
sont convergentes pour toutes les valeurs de la variable ¢ et que les intégrales © 
dans les derniéres lignes des équations (13a) et (13b) tendent vers zéro lorsqu’on 
pousse l’itération a l’infini. Pour cela, il n’y a qu’a considérer le systéme d’équa- 
tions majorantes 


ad? z d*y 
a =P (lts+ 9) =X, pap tee yee eae 


ot p est le plus grand des nombres Bl, er 
En suivant la méthode du numéro précédent, on a d’abord 


ae nef Ot at yan set aty—p ft ( Ge + GP) a 


BO epg ary Eka dy Pup : 
=1X ee (a tar)+ 5 fex+ ya, 
dott 


/ »\ d » dy : 
(1 zs F 1) > 4 Le ie + * i) .e* £ 8 Y 


Or F fe (S Se de, 


En substituant aux dérivées de X et de Y leurs expressions au moyen des — 
dérivées de x et de y, tirées de (14), ona 


A eed DL may 4 C3 
+50) + ee g-(+fajx+Feay 
2 
7 Se 0 (Se + ) dt 


et, d’une maniére générale, le procédé d’itération donne 
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(1 + P t2 4 2p? pee ! 2 p> AD ese S00. Hl ae are Pe pen) 


2! 4) 6! (2 n)! at 
(ee a Sn Oe 

= (« ek -- 4 eee go fovee 4 Ow i a 7 x 
+ (47 oe 5 Si a gear e 9) ee = fe ( XY eae 


On a ici, l’intégrale étant prise entre des limites finies, 


.  gn-l pn 
lim 27)! 


nC 


2m (X 4+ VY) dt=0. 


D’autre part, la limite du rapport des termes de rang (m + 1) aux termes 
de rang ” tendant vers zéro avec 1/n, quelque soit ¢, les séries qui figurent comme 
coefficients de dx/dt, dy/dt, X et Y sont convergentes et d’une convergence 
rapide pour toutes les valeurs de ¢. A fortiori les coefficients correspondants 
dans les formules (13a) et (13b) convergent pour toutes les valeurs de ¢, et on a 


eee Sky, Bo CD ee X ey YS (15) 
avec 
A=1+S 8+ pie , a=tt ay fhe 9 
Be par NE dee sat ~~ pe spits aay nes 
egy ret Sart 
Les équations (15) donnent enfin 

fae Sigh blak =p = lak tle (17) 

avec 
ee ee ht ee 


Observons d’abord que l’on a 


AC— Bta=1 +7384 


48 KyriILLE Poporr ZAMP 


Pour montrer qu’on a Ly, = Ly, il n’y a qu’a former la différence 
(CB — By) — (AB — Ba) = B(C — A) — Bly—a@). ; 


Un calcul simple, en tenant compte de (16), montre qu’elle est nulle. \ 

Dans ce qui précéde, nous avons mis nulles les constantes d’intégration, 
ce qui revient 4 admettre qu’aux forces nulles correspondent des flux nuls, 
comme c’est admis dans les différentes théories des phénoménes physico- 
chimiques. Dans le cas général, il faut écrire 


AL BD aaX+PY+M, B14 CD apX+yV4N, 
} 
M et N étant des constantes d’intégration. . 
Remarque. On arrivera aux mémes résultats en prenant comme point de 
départ une autre fonction caractéristique au lieu de l’entropie, ce qui généralise_ 
les résultats concernant les transformations linéaires des forces et des flux 
étudiées dans l’ouvrage cité de DE GROOT. : 


Expression des coefficients L;; au moyen de a, b, c ; 


Les relations phénoménologiques 


o* = A(ax+by) + B(bxtey) =. Be SA 


(A, B, C des constantes) 
P= Blaxt+by)+C(bxt+cy) =BX+CY, | 


ary 
\O 
— 
ee 


ou, pour éviter les indices, nous avons mis A, B, C au lieu de Ly,, Ly, = Loy, 


Ly. respectivement, conduisent, par un choix convenable de A, B, C, aux équa- 
tions 


d?x : d*y 
ae a = ve oe 
En effet, des équations (19), on tire immédiatement 
13, 
aie ~ (40+ Bb) 94 (4b4 BH) @ 


=[(Aa+Bb)A+(Ab+ Bc) BX 
+[(Aa+ Bb) B+(4b+ Bo) C]=4,X+B,Y, 


=([(Ba+Cb)A+(Bb+Cc) BX 


| 
| (21) 
J 


+ [(Ba+Cb)B+(Bb+CoC]}Y=D,X+C,Y. 
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On constate immédiatement qu’on a D,=B,. On démontre ensuite 
qu’on peut déterminer A, B, C comme fonctions de a, b, cde facon qu’on ait 
/D,=B8,=0 et A4,=C,=pu=1. En effet, en écrivant B, =0 et Det) 
sous la forme 


Aat+Bb=—-<(4b+Bo), Bee == (Ba Cb 


by] ds 


Péquation A, = C, = uw donne 


Avo + Bie 


B (B2— Ac) 3B (BSAC i= 1 
d’ot 
Ab+Bc=Ba+Cb; C=A+B—“. 
_ Avec cette valeur de C, l’équation A, — C, = 0 devient 
A \2 A\ ¢ CG 1a) 
ee a ele = 


d’ot 


A c ac— b? C a Geo 
poo |) Le oe =e pre? 
ce qui donne enfin 


p= sr (we— DB) [a+cF2Vac—b]. 


En choisissant Bde faconqu’onait « = 1, les équations (21) deviennent enfin 


it ace eal 2) eee 


ce qu'il fallait démontrer. 


Les relations phénoménologiques (19) représentent un systeme 
d’intégrales premieres du systéme (20) 


Ecrivons les équations (20) sous la forme 


d ’ dé dy dn 
ee ene OOD ate OY ag oe ee. 


Elles admettent les intégrales 


: an t 
Rae ene. € =H 6 es Heer, 
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ou y est déterminé par l’équation 
v4—r2(a+c)t+ac—b?=0 


et «, 6, y sont donnés par 
y=ra=Pa. 


p=r, “g= — oe, 


Les racines de ]’équation en ¢ sont ici 
c Ce Cla sae 
+ VE) +b%?—ac , h%=—nN, 


i 
At eg ters 


ede eed ee 
La forme quadratique a x?+ 2b x y+ c y* étant définie positive, on a. 
fick QO eee 


b> —ac= 0, 


et par conséquent 


On a ainsi 
FP gee Ue Oe. ay, Fo 


Les intégrales du systéme (22), s’annulant pour ¢ = +-0e, sont par conséquent 
ye Coe + Cee. . 
(23) 


RCo eth Owe, 
——— rot ryt ‘° 
f= Cofre+ CBee, n=Cyyge™ + Cyyge™. 


En substituant a e”', e”' dans les derniéres équations leurs valeurs tirées 
des deux premicres, ayant égard aux relations (23), on obtient aprés des calculs 


y=A(ax+by)+ B(bx+cy) 


simples 
ax» _ Ba%— Bia, | By — By 
= = x — 
dt 1, — Oy h, — Oy 
=AX+BY=(Aa+Bb)x+(Ab+Bo)y, 
seas a Me °, te ioe 
ge me Daerah) othe Colby aoa 


dy _ a 
tlie. Xe —'K, 


dt 
=DX+CY=(Da+Cb)x+(Db+Co)y. 
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La comparaison des coefficients de x et y dans les premiers et les seconds 
membres de ces équations, en égard aux expressions de 7,, 74, «;, B;, y;, donne 
apres quelques transformations 


we 6) B=" (7,=7,) Vac~— b2, (ac 6*)D==(r,= 4) Vac— 67", 


dowB = D, (L,.= L,,), ce qu'il fallait démontrer. 


Summary 


Basing on the probability theory, ONSAGER and Casrimrr have established 
that the coefficients L,;, of the phenomenological relations between the ‘‘ Flows”’ 
J; and the “Forces” X;, J; = 3’ Lj, X, near the steady states [i.e. taking the 

z 


expression of variation of entropy as AS = —()) g;, ¥; ¥;,)/2] satisfy the equalities 
Ly, = as Lig. cee 

' In the preceding it is shown that these are phenomenological consequences. 

of RESTRON!S Mechanics in the more general case where: AS is ‘assumed as 


ous Bj Hy = (x Sin ¥4 Xp) /2. 
t 1, 


(Recu le 14 juillet 1951.) 


Lagrange-Hermitesche Interpolation im Komplexen’) 


Von Hetnz UNGER, Darmstadt?) 


1. Einleitung 


Die Interpolation einer analytischen Funktion im Komplexen wird mit zu- 
nehmender Bedeutung nichtelementarer Funktionen ein immer dringlicheres 
Problem. Die dabei zu beachtenden Gesichtspunkte unterscheiden sich im all- 
gemeinen wesentlich von denen der reellen Interpolation. Die Tafeln erfordern 
im Komplexen doppelten Eingang und kénnen kaum so engmaschig angelegt 
werden, da etwa kubische Interpolation [1]*) noch ausreicht. Es sind von 
vornherein hochgradige Interpolationsformeln ins Auge zu fassen. Dabei muB 
man versuchen, mit mdéglichst wenig Angaben auszukommen, um die Tafel- 
werke nicht zu umfangreich werden zu lassen. Andererseits mu8 aber auch der 
Aufwand zur Berechnung eines Zwischenwertes mit hochgradigen Naherungs- 


1) Es handelt sich um eine ausftihrliche Wiedergabe eines Vortrages, der aut der Tagung der 
Gesellschaft fiir angewandte Mathematik und Mechanik (Freiburg i. Br. vom 29. bis 31. Marz 1951) 
gehalten wurde. Vgl. auch Z. angew. Math. Mech. 31, 246-247 (1951). 

2) Institut fiir Praktische Mathematik der Technischen Hochschule Darmstadt. 

3) Die Ziffern in eckigen Klammern beziehen sich auf das Literaturverzeichnis auf S. 64. 
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polynomen in ertraglichen Grenzen bleiben. Beide Forderungen richtig gegen-. 
einander abzuwagen, rechtfertigt besondere Untersuchungen. Es handelt sich 
hier vornehmlich um die Berechnung einzelner Zwischenwerte, nicht um die 
Aufgabe der Untertafelung, die wieder andere Methoden erfordert. 

Ineiner Arbeit von SCHNELL [2] ist gezeigt worden, daB bei der Interpolation 
im Komplexen die Aufspaltung der Funktion in Real- und Imaginarteil ein 
Umweg ist. Die Verwendung der bekannten Interpolationsformeln, wie Bessel- 
sche, Stirlingsche, Everettsche Formel, mit komplexem Argument kommt 
gleichfalls kaum in Frage, da die Stiitzwerte bei hochgradiger Interpolation 
im allgemeinen zu weit vom interessierenden Bereich abliegen, was eine un- 
giinstige Restgliedabschatzung zur Folge hat. 

Als naheliegende Stiitzstellenanordnung erscheint im Komplexen das regel- 
maBige Polygon. Dabei gelangt man zu besonders vorteilhaften Formeln, wenn 
man nicht, wie bei der itiblichen Lagrangeschen Interpolation, allein die Funk- 
tionswerte in den » Eckpunkten des Polygons heranzieht, sondern nach der 
auf HERMITE [3] zuriickgehenden Erweiterung der Interpolationsformeln auch 
noch die Ableitungen bis zu einer bestimmten Ordnung beniitzt. Man kommt 
auf diese Weise schon mit wenig Stiitzstellen zu Interpolationsformeln hohen 
Grades, wie sie fiir unsere Zwecke bendtigt werden. Wir setzen fiir das Folgende 
voraus, daB in den Polygoneckpunkten die Funktionswerte und die g ersten 
Ableitungen (in allen Punkten bis zur gleichen Ordnung?*)) bekannt sind. Die 
Polygoneckpunkte bilden dann die «Gitterpunkte», fiir deren Argumente die 
Funktionswerte nebst Ableitungen zu tabellieren sind. 

Es wird sich zeigen, da man sich dabei unter Verwendung eines kartesi- 
schen Koordinatensystems im allgemeinen auf das Quadrat als besonders 
zweckmaBigen Sonderfall des regelmaBigen Polygons beschranken kann. 

Fir die praktische Durchfiihrung der Interpolation miissen die Interpola- 
tionskoeffizienten in einer Hilfstafel vorliegen. Im Falle des Quadrats aber 
brauchen von den (q+ 1) bendtigten Koeffizienten nur g+ 1 in einem Quadrat 
der Seitenlange 1 mit einer bestimmten Maschenweite (zum Beispiel 1/100) 
vertafelt zu werden. Den gleichen Vorteil kann man bei anderen Polygonen 
(zum Beispiel Sechseck) nur unter Verwendung von Polarkoordinaten erreichen. 
Im allgemeinen ist aber das Arbeiten mit kartesischen Koordinaten bequemer. | 

Die Formeln werden zunichst alle so aufgestellt, daB der Polygonmittel- 
punkt als Stiitzstelle unberiicksichtigt bleibt, da dies den haufigeren Fall der 
Stiitzstellenanordnung darstellt. In dem erganzend behandelten Fall mit Be- 
riicksichtigung des Mittelpunktes miissen dann auBer den g +1 Koeffizienten 
der Eckpunkte noch die Koeffizienten fiir den Mittelpunkt bereitgestellt wer- 
den; doch braucht dies nur in einem Teile des Quadrats zu geschehen. 


5 Der allgemeinere B all, daB OAT Ableitungen IM 21, gy in 2) usw. berticksichtigt werden, wird 
hier nicht betrachtet. Weiter bleibt die Verallgemeinerung, daB Ableitungen von beliebigem Grade 
herangezogen werden, einer besonderen Arbeit vorbehalten. 
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Den Vorteil der Kenntnis nur einer kleinen Zahl von Interpolationskoeffi- 
zienten kann man bei regelmaBigen Polygonen auch dadurch ‘erreichen, daB 
man Zum Beispiel auf der Mitte jeder Polygonseite eine weitere Stiitzstelle hin- 
zunimmt und mit der gewohnlichen Lagrange-Formel arbeitet. Es erscheint 
aber im allgemeinen wesentlich einfacher, an den Eckpunkten die Ableitungen 
einer speziellen Funktion mit heranzuziehen, da diese oft durch Rekursions- 
formeln oder Differentialgleichungen gebildet werden k6énnen. SchlieBlich ist 
auch der Aufbau der Interpolationskoeffizienten nach LAGRANGE-HERMITE 
besonders einfach, was fiir die Berechnung dieser Koeffizienten eine nicht un- 
wesentliche Rolle spielt. 

Vorliegende Arbeit beschaftigt sich mit dem Aufbau der Interpolations- 
koeffizienten und ihrem Zusammenhang untereinander. Beziiglich der prak- 
tischen Durchfihrung sei auf [4] verwiesen. 


2. Lagrange-Hermitesche Interpolation 


Wir setzen voraus, daB an den durchweg als verschieden angenommenen 
Stiitzstellen 


Cah Lay erate wei (1) 


der komplexen Veranderlichen z = x + 7 y die Funktionswerte und die g ersten 
Ableitungen der zu interpolierenden analytischen Funktion g(z) gegeben sind: 


FE Ege snd fees ane Sem ees operas tan (2) 

Damit liegen  (¢g +1) Angaben vor, und wir suchen ein Polynom P,(z) vom 
Grade y mit 

ye og 1 1) (3) 


welches an den Stiitzstellen die vorgegebenen Funktionswerte und Ableitungen 
annimmt. 


Wir setzen 
m 4 
Dae ALE 2) Ge (4) 
p=1p=0 
und damit 
AO TAG cme) (5) 


mit dem Restglied R,,4(z). Zur Vereinfachung der Schreibweise lassen wir im 
folgenden bei den Interpolationskoeffizienten Ea eS die oberen Indizes n, q 
wieder fort, da Verwechslungen kaum zu befiirchten sind, schreiben also ein- 
fach L,,p(z). Wir filhren folgende Vektoren ein?): 


1) Die geschweifte Klammer wird zur bequemeren Schreibweise einer einspaltigen Matrix 


verwendet. 
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: , ie : (g) < (q) 
g ={ 8i,ds ewe Bixee bnr ares pe eX! 
pr Re Lend se 


34 = eee, We 
Weiter beniitzen wir die [7 (¢ + 1)]-reihige Matrix 
a3, Ap a3, “| (7) 


ea 1 2 e's 

Brg = [ar dar 38n day ot Ciedage tea, wae wake 

deren Determinante eine Verallgemeinerung der bekannten Vandermondeschen 

Determinante det %,, darstellt [5] und sich durch diese wie folgt ausdriicken 
1aBt?) : 


q(q+1) n(n —1) 


det B,g=(-D 2 SLI gl aidet ial (8) 
mit 
1,n 
det Ba = [] (. — 4). (9) 
por : 


Sind also, wie vorausgesetzt, die z, voneinander verschieden, so muf mit 
det Y,, auch det B,,., von Null verschieden sein. 
Setzen wir daher fiir das Naherungspolynom an: 


Pale; c= - ay es (10) 


so erhalten wir zur Bestimmung der y + 1 Koeffizienten a, das Gleichungs- 
system: 


Base C=; (11) 


welches eine eindeutige Lésung besitzt. 


In Erweiterung der Darstellung der Lagrangeschen Interpolation [6] laBt 
sich P,(z) in der Form schreiben: 


2 LP G 
P,(z) = — det | f q 4 nya cae | ea (12) 
wahrend das Restglied lautet: 
) B, € 
Ne 1(2) == det ( 3’ @ a det Pe os (13) 


1) Auf die Literaturstelle [5] wurde ich von Herrn Prof. Dr. U. WEGNER in dankenswerter 
Weise hingewiesen. Dort wird ein rekurrentes Verfahren zur Berechnung der verallgemeinerten 
Vandermondeschen Determinante angegeben, die sich von (7) nur durch andere Zeilenanordnung 
und unwesentliche Faktoren unterscheidet. Der allgemeine explizite Ausdruck findet sich dort 
nicht. Das Ergebnis (8) wurde von mir auf andere Weise gefunden, indem eine gewohnliche Vander- 
mondesche Determinante von der Ordnung n (q + 1) differenziert wird. Auf diese Weise laBt sich 
auch der allgemeine Fall, vgl. FuBnote S. 52, mit q, Ableitungen in ay leicht erledigen. 
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Zur Durchfithrung der Interpolation werden die 2 (g + 1) Polynome L,, ,(2) 
vom Grade y 
¢ L  5(2) as 3 Cup ; (14) 
mit 


Cup ={(G,,s)0. (C, p)a» prac: (C,»)y} (15) 


benotigt. Zur Bestimmung der Koeffizienten (Cu,p)» Sind 1 (¢ +1) Gleichungen 
erforderlich, die mit dem Vektor 


e519 Oe raise L AO ate st 4 (16) 


v-te Zeile 
folgendermaBen lauten: 
Wile Cp Ones ie ( (=e 1 2iee 8 pi Ol. O\iA (12) 


_ Die Bestimmung der Koeffizienten ist unter den gemachten Voraussetzungen 
eindeutig. Fassen wir die » (g + 1) Vektoren c¢,,, zur (y + 1)-reihigen quadrati- 
schen Matrix 


tC od Corus on Cain etebicge Sees Gag) (18) 


zusammen, so folgt aus (17): 


Bi Oa. a = €. (19) 


Die Matrix der Koeffizienten ist also die Kehrmatrix?) zur verallgemeinerten 
Vandermondeschen Matrix (7) (vgl. auch [7]). ° 

Da sich die Vandermondesche Determinante (9) nicht andert, wenn z durch 
z — z* ersetzt wird, so dndert sich auch det %,,, nicht. In den Spalten der 
Matrix ©,,,, deren v-te Zeile die (vy —1)-te Ableitung von L, ,(z), dividiert 
durch (vy — 1)!, enthalt, tritt dieser Ausdruck dann an der Stelle z= z* an- 
statt in z= 0 auf. 

Fiir die zu interpolierende Funktion g(z) miissen wir voraussetzen, daB g(z) 
im Innern und auf der Berandung des kleinsten konvexen Polygons, das durch 
die Punkte z, 2,, 2), ---, 2, bestimmt wird, regular ist. Allerdings muB fiir die 
praktische Durchfiihrung zusatzlich gefordert werden, daB die nachste Sin- 
gularitat noch geniigend weit entfernt liegt, um eine brauchbare Restglied- 
abschatzung zu erhalten. 

Von dem Restglied, das in allgemeiner Form wie folgt lautet: 

1 ¢ dt 

BZ) = [2 — 2)... (2 — ae zat 1 oe ae ae —2z,)]!  f—2 


, (20) 


1) Einen Spezialfall dieser Beziehung mit einer etwas abgeanderten verallgemeinerten Vander- 
mondeschen Determinante findet man bei KowaLrwsk1 [7]..Die Elemente der Kehrmatrix stellen 
auch die Koeffizienten der Partialbruchzerlegungen gewisser leicht angebbarer rationaler Funk- 
tionen dar, worauf schon in [5] hingewiesen wird. 
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wobei C den Rand des oben erwahnten Polygons darstellt, interessieren uns_ 
in diesem Zusammenhange hauptsachlich Abschatzungen. Dazu wird die Formel 
(20) unter Anwendung bekannter Beziehungen (vgl. [6], S. 16, Darstellung des 
Restgliedes nach HERMITE, Darbouxscher Mittelwertsatz usw.) umgeformt. 
Hierbei sei auch auf die Arbeit von P. MonTEL [8] verwiesen. 


3. Einheitspolygon in der komplexen Ebene 


Wir legen als Stiitzstellenanordnung ein regelmaBiges Polygon von m Ecken 
und einer Seitenlange 1 zugrunde, im folgenden kurz Einheitspolygon genannt, 
und wahlen ein z-Koordinatensystem mit dem Nullpunkt im Polygonmittel- 
punkt. Die Numerierung der Eckpunkte z, sei so durchgefiihrt, daB 


ih a 
1 as G \ 4 
= = (21) 
aaa 2 sin — 
é — Il r 
und daB 
t —~ (4-1) 
2, = %e F (22) 
so daB 


wird (Figur 1). Die Wahl dieses Polygons bedeutet keine Einschrankung, da 


2 


Fig. 1 


Lage des Einheitspolygons. 


jedes regelmaBige Polygon sich auf das gewiinschte durch eine Lineartrans- 
formation zuriickfithren 1aBt. Laut Voraussetzung soll sich die zu interpolie- 
rende Funktion g(z) in und auf dem Rande des Polygons regular verhalten. 
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Die Vandermondesche Determinante det %,, hat den Wert: 
det G,, = (2sin =)" a ATE iz teat) sal 


_ Bei der Quadratwurzel ist immer das positive Zeichen zu nehmen. 
| Neben der Veranderlichen z werden wir noch die Veranderliche 


W=2— 2, (24) 


_ beniitzen, deren Achsensystem in Figur 1 gestrichelt eingetragen ist. In ihr 
sind die Koordinaten der laufenden Eckpunkte des Polygons dann: 


—1 : 
(Q), = URS Hg I, on (25) 


Die Veranderliche z bzw. w beschranke sich auf das Innere und den Rand 
des Einheitspolygons. 

Das Einheitspolygon stellt die fiir die komplexe Ebene natiirliche Verall- 
gemeinerung gleichabstandiger Stiitzstellen auf der reellen Achse dar. 


4. Die Lagrange-Hermiteschen Koeffizienten fir das Einheitspolygon 


Setzt man die Interpolationsformel (5) fiir das Einheitspolygon an, so muB 
der zu berechnende Wert unveradndert bleiben, wenn man eine Lineartransfor- 
mation von z durchfiihrt. Macht man insbesondere die Drehung 


2% 


2) zen" tay NS 28S testy (26) 


bei der der Punkt z= z, im Uhrzeigersinn in den Punkt z, tbergefihrt wird 
(vgl. Figur 1), so zeigt sich, daB von den m (¢+1) bendtigten Koeffizienten 
nur g +1 bereitzustellen sind. Denn es gilt nach Einsetzung von (26) in (5) 
die Beziehung: 


De net cw eae (aye cen) (27) 


Fiir g = 0 erhalt man fiir die gewohnlichen Lagrange-Koeffizienten, die wir 
mit L,,(z) bezeichnen, den Zusammenhang: 


L (2) = 2 Bal 


Fiir jeden p-Wert wird also nur die Kenntnis eines Koeffizienten, namlich 
L,,p, im Einheitspolygon an m Argumentstellen bendtigt. Andererseits wiirde 
es auch geniigen, jeden Koeffizienten in dem m-ten Teil des Polygons zu kennen. 
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Interpoliert man im Reellen mit aquidistanten Stiitzstellen, so bendtigt man, 
wenn 2 gerade ist, (¢+1) /2, bei ungeradem n dagegen (g+1) [(2+1)/2] 
Koeffizienten. Hierauf ist fiir die gewdhnlichen Lagrange-Koeffizienten zum 
Beispiel in [9], S. XVII, hingewiesen. Im Einheitspolygon werden die Verbaleg 
nisse also besonders einfach. 

Ubrigens erhalt man durch Differenzieren von (27) fiir die k-te Ableitung?) 


(u—1) (p—k) 


. 22 
EY (2) = Lee 


Hs; p) 
womit auch die Ableitungen auf diejenigen mit dem Index yw = 1 zuriickge- 
fiihrt sind. 

Die Interpolationsformel (5) lautet unter Beriicksichtigung von (27): 


yo see 
gz) = > SY Ly glee» PY ol) 4 R (2), (5a) 


wobei 2”) nach (26) zu bilden ist. ee. 

Fihren wir an Stelle von z die Veranderliche w nach (24) ein und bezeichnen 
wir die Lagrange-Hermiteschen Koeffizienten als Funktionen von w mit M u,p(@) 
und die zu interpolierende Funktion mit G(w), so geht (5a) iiber in 


ms q , ow 
= 37S My sal et OP OP 4 RE, 0) at 
u=1 p=0 
mit 


2% 
i (1~ 4) a3 
w) —we n + We gs SUE eae a ee (28) 


n—f i 
Fur die praktische Durchfithrung der Interpolation wird im allgemeinen die 
Veranderliche w bevorzugt (vgl. [1]). 
Aus (5a) folgt die fiir Kontrollen wichtige Beziehung 


ok n k . 27 
BO" Ss (“1 >! p (u—1) 
ay =D, dy Tole he 

P 


B=1 0 


die fiir die gewohnlichen Lagrange-Koeffizienten mit k = 0 iibergeht in 


" 


Le)” Life). 


Mel 


Da die Lagrange-Hermiteschen Koeffizienten L,,,p(2) an allen Stiitzstellen 
auBer in z= 2, eine (g + 1)-fache Nullstelle peslteen miissen, laBt sich ein | 
wesentlicher Bestandteil durch die gew6hnlichen Lagrange-Koeffizienten aus-— 


1) Unter Li om ist jedesmal die Ableitung nach dem Argument zu verstehen. 
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-driicken. Es gilt die Beziehung: 


Ly, p(2) = [E,(2)7"1 4, 52). (29) 


Hierbei stellt /,,,(z) ein Polynom g-ten Grades dar. Von den n (¢+1) Glei- 
chungen (17) sind durch [L,,(z)|**1 bereits (w — 1) (q¢ + 1) erfiillt. Die restlichen 
Gleichungen 


(k) ee (ure : 
Veet Gris be Bh it 051523, cog (30) 


-dienen zur Bestimmung der qg +1 Konstanten des Polynoms J, ,(z). Dieses 
Polynom mu8 daher noch durch (z — z,)?/p! teilbar sein, woraus der fiir # = ¢ 

_ besonders'einfache Aufbau der Koeffizienten folgt. 

Die Beziehung (29) gilt fiir jede Stiitzstellenanordnung. Beim Einheits- 
polygon wird sie infolge (27) nur fiir «4 =.1 bendtigt. 

Eine Restgliedabschatzung ist bei der komplexen Interpolation und ins- 
besondere bei der hier beniitzten hochgradigen von noch gréBerer Bedeutung 
als im Reellen, da hier die Genauigkeit eines Zwischenwertes nur durch eine 
sorgfaltige Restgliedabschatzung garantiert werden kann. 

Bedeutet K,j(¢,3) eine obere Schranke fiir den Betrag der [ (¢ + 1)]-ten 
Ableitung der Funktion g(z) auf dem Rande des Polygons, so 1aBt sich folgende 
Abschatzung des Restgliedes angeben: 


ene n “q+ 
(cos =) tee : Kn (q+1) 
n(q 
| Ree a2) | < " (2sin2)" in (q+ 1]! ° (31) 
Wainer 


Die Hauptaufgabe bei der praktischen Fehlerabschatzung ist die Abschatzung 
der [m (q + 1)]-ten Ableitung auf dem Rande. Fiir eine Orientierung geniigt es 
oft, die Abschatzungen in den Eckpunkten zu kennen. 


5. Mitnahme des Polygonmittelpunktes 


Nimmt man den Mittelpunkt des Polygons als weitere Stiitzstelle hinzu und 
setzt dort wieder den Funktionswert g\?) als bekannt voraus, so erhalt man an 
Stelle von (5) folgende Interpolationsformel : 


n q q 
gz) = D7 DN ol?) ef) + 2, Noval) go” + Ryesal2), (5c) 
w=1 p=0 p=0 


-wobei y* = (n+1) (¢+1)—1 den Grad der Interpolationskoeffizienten dar- 
stellt!), Fithren wir wieder die Transformation (26) durch, so erhalten wir fiir 


1) Wir lassen auch hier die beiden Indizes 1 -- 1 und q weg und bezeichnen die Koeffizienten 
mit Ny,p(z), weil keine Verwechslungen zu befiirchten sind. 
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die Polynome N,,,(z) die zu (27) entsprechende Beziehung, wobei also fiir 
jeden #-Wert ein Polynom benotigt wird. AuBerdem gilt: 
Ny ol@) = (3) Ex) By (2) = ENG) U8) 
iP : 
Dabei sind auch die Polynome /#,(z) wie J;,,(z) vom Grade g, und N,(z) ist der 
gewohnliche Lagrange-Koeffizient. 
Fir die beim Mittelpunkt auftretenden Polynome gilt: 


2x 
i—p(u-1) ., 
No, p(2) = No,,[2]e * fee far = 2,57 


Diese Koeffizienten brauchen also nur in einem Dreieck (1/n des Polygons) 
bekannt zu sein. Unter Heranziehung des gewohnlichen Lagrange-Koeffizien- 
ten N4(z) kann der Koeffizient wie folgt ausgedriickt werden: 


a ap : in (n/2 ere SE oe 
No, p(2) = [No(z)]2*1 pr mit N,(z) = e™) (2 sin =) 2e+ 1, 


Bei Vorliegen einer Potenztafel im Komplexen ist dieser Koeffizient einfach 
zu bilden. 

Bei Mitnahme des Mittelpunktes empfiehlt es sich, die Koordinate z auch 
fiir die praktische Rechnung beizubehalten. 

Es sei noch darauf hingewiesen, daB auch die Formeln, bei denen das 
Polygon um z/n gedreht ist, sowie der Fall, daB die Stiitzstellen mit Funktions- 
werten und Ableitungen auf der Mitte der Polygonseite liegen, leicht auf dem 
hier beschriebenen Wege hergeleitet werden kénnen. 

Fir das Restglied gilt mit entsprechenden Bezeichnungen wie bei (31) die 


Abschatzung: 
| hae i ma \n “Tq+1 
(cos =) + 1 x 


(ne) (+0) 


tes 7 (2 sin =)" [(~ + 1) (g+1)}! 


6. Das Quadrat als Einheitspolygon (ohne Mittelpunkt) 


é 


Fiir die praktische Durchfiihrung der Interpolation erweist sich die quadra- 
tische Stiitzstellenanordnung als besonders vorteilhaft. Hier hat man die Koef- 
fizienten L,,,(z) fiir p= 0, 1,..., q zu z-Werten eines quadratischen Netzes im 
Einheitsquadrat, also einem Quadrat mit der Seitenlinge 1, zu vertafeln. 
Wahlt man als Maschenweite des Netzes zum Beispiel 1/100, so sind fiir jeden 
der g +1 Interpolationskoeffizienten L,,y(2) 10000 Werte anzugeben. (Durch 
Spiegelung an der Diagonale von 2, nach 2; kénnte diese Zahl noch auf etwa 
die Halfte reduziert werden, vgl. hierzu [4].) Ist z dann ein Tafelargument, so 
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sind beim Quadrat auch die Werte z nach (26) fiir 4 = 2, 3, 4 wieder Tafel- 
argumente, die zu z in ganz einfacher Beziehung stehen, namlich: 


20 —xtiy, 2<2ay—ix, 2) =2—x¢-—1y, 222 —ytin. 


Damit kénnen-die Koeffizienten L,,, fiir alle vier Argumentwerte 2 bequem 
aus der Tafel abgelesen werden. Die Berechnung der L,,,(z) erfordert nach (27) 


Pe ee ee ice 2,334 (27a) 


dann nur noch die Multiplikation mit Potenzen von 7. 

Der Ubergang auf die fiir die Vertafelung bequemere Veranderliche w be- 
reitet keine Schwierigkeit. 

Fiir ¢g =.0, 1 und 2 ist die Vertafelung der Koeffizienten L,,, in der Ver- 
anderlichen w zur Zeit im Gange, und zwar mit einer Maschenweite von 1/100, 
die sich aus Restgliedabschatzungen fiir die meisten Zwecke als die giinstigste 
ergab [4]. 

Will man feiner als mit der in der Interpolationstafel vorgegebenen Maschen- 
weite (zum Beispiel 1/100) unterteilen, so kann man entweder zundachst die 
vier umliegenden in der Tafel enthaltenen Eckpunkte in der angegebenen 
Weise berechnen und dann mit einer Interpolationsformel niederen Grades 
nachinterpolieren, deren Genauigkeit dann stets ausreichen wird; oder aber 
die Interpolation wird nach anderen Gesichtspunkten durchgefiihrt, worauf 
teils in Abschnitt 7, teils in einer spateren Arbeit eingegangen wird. 

Zur Berechnung der Funktionswerte in der Quadratmitte erhalt man fiir 
q = 0, 1 und 2 die folgenden Formeln: 


0) = = [g. + £2 + bs + Bal + Ra, | 


1 
4 
il - 5 i y ° / / ° / 
Be 7 eit faa est Salt ep Cle th = byt eal Ky, 
(32) 
1 - 5 + IZ ° [2 / ie / 
lita ba bs Calte —png 180 282 6s 18a) | 
+ see (-e tes —el+Re. | 
Dabei sind die Funktionswerte und Ableitungen an den Stellen 


14% is 1+2 Pew 
a> ee Fie Doe 2 


B= 


zu nehmen. 
Aus (31) erhalt man hierbei fiir Rj, die Abschatzung 


Aerie ni Oren 


LJ 
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Dies bedeutet, daB zum Beispiel die Zylinderfunktion /,(¢) fiir €¢=5,5+ 1 5 
nach (32) aus den Werten fiir : 
: 


= 544, Co=64+4, Cg=6421, y= S+2t 


auf neun Dezimalen genau berechnet werden kann, wenn die Stiitzwerte sowie 
die ersten und zweiten Ableitungen auf neun Dezimalen genau bekannt sind. 

Vergleicht man bei quadratischer Stiitzstellenanordnung die Interpolation 
nach LAGRANGE-HERMITE unter Beniitzung der ersten und zweiten Ablei- 
tungen mit der Interpolation nach der gewohnlichen Lagrange-Formel, so wird 
letztere im Mittel bei gleicher Genauigkeit etwa 1/100 der Maschenweite er~ 
fordern, angenahert gleiche Betrage fiir die im Restgled auftretenden Ablei- 
tungen vorausgesetzt. Man spart also durch die Verwendung der Lagrange-_ 
Hermiteschen Formel in einem einzigen Quadrat die Angabe von etwa 10000 
Funktionswerten bei einem Aufwand von 12 komplexen Multiplikationen zur 
Berechnung eines Zwischenwertes, der in Anbetracht der groBen Ersparnis an_ 
Zahlenmaterial wohl vertretbar ist. : 

Auf besondere Eigenschaften der Koeffizienten L,,(z) wird in [4] einge- 
gangen, wo auch besondere Kontrollformeln aufgefiihrt sind. 


7. Eine Umschreibung der Lagrange-Hermiteschen Formel 


Zu einer anderen Vorgehensweise bei der Auswertung der Lagrange-Hermite- 
schen Interpolationsformel — hierbei erweist es sich als zweckmaBig mit der 
Variablen w nach (24) zu arbeiten — gelangt man, wenn man bei der Formel 
(5b) aus dem Koeffizienten M,,,(w) den von p unabhangigen Bestandteil 
[M,(w)|** gemaB der Beziehung (29) herauszieht: 


M,,,(w) = [M,(w)]*** A, ,(w), (29a). 


worin /,,,(«) wie das entsprechende /, ,(z) ein Polynom g-ten Grades ist: 


G(w) SIM, [w') eye 5 Ay p[w(] e'? nn “-) GO) 4 Re (w). (Sb*)N 


7 
u=1 


Die Koeffizienten von A,,,,(w) ergeben sich aus der fiir die Verinderliche we 
umgeschriebenen Gleichung (30) und kénnen fiir ein bestimmtes und g ein 
fiir allemal berechnet werden. Wir setzen 
‘ 1 
Ay, p(w’) a 2) a), p w" mit ay,» = p! (4,, 0 = Oatorsy = p) (33) - 


und fassen die a, , in einer (¢ + 1)-reihigen quadratischen Matrix zusammen: 
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1 0 Claas me Fert), 
1 
ayy Ty 0 0 
1 
Wr, CS 0 ; (34) 
1 


q! 


feet 0 Cw 
9 ; 
ee = (1-3) i 0 
2 9 F A 
a 


Ordnet man die zweite Summe in (5b*) um nach Potenzen von w bzw. w'), so 
erhalten wir: 


n qd 
G(w) =] om [M, [wo] q+1 oy. a, [wo -- i), 4 (w) K (5 b**) 
ual v=0O 
Die Koeffizienten d,, , ergeben sich aus: 
; dig Ex9 Gy 
(1) 
dix LaDs| Gi, : i v(ye—1) 
; : MMe tej” 
unde ve — ON e.g 
9) 
dig Eng G, 


Im Falle des Quadrats gehen die Faktoren e,,,, wieder in Potenzen von 2 tiber. 
Die Dreiecksform der Matrix Y,,,, vereinfacht die Ausrechnung der Koeffi- 
zienten d,, , (bei Verwendung der Veranderlichen z statt w hatte man eine volle 


Matrix erhalten). Insbesondere gilt: 


Die Koeffizienten d,,, enthalten die Funktionswerte und die q ersten Ab- 
leitungen. Man hat also fiir eine zu interpolierende spezielle Funktion G(w) 
an den Eckpunkten des Einheitspolygons diese Koeffizienten bereitzustellen. 
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Will man aus diesen Werten d,,, umgekehrt die Ableitungen G{?) selbst wieder! 
errechnen, so kann dies infolge der Dreiecksgestalt der Matrix U,,,z leicht durch: 
Aufrechnung des gestaffelten Gleichungssystems erfolgen. 

Liegt die spezielle Funktion G(w) in dieser Weise vertafelt vor, so wird die: 
Interpolation nach (5b**) durchgefiihrt. Dabei ist es zweckmaBig, die zweite: 
Summe bei Nichtvorhandensein einer Potenztafel mit dem Horner-Schema aus- : 
zuwerten. Als einziger Koeffizient miiBte dann noch [M,(w)]*** im Einheits-. 
quadrat vertafelt vorliegen (bei der friiher gebrauchten Vorgehensweise muBten | 
g +1 Koeffizienten vertafelt werden). Da aber dieser Koeffizient auch verhalt-. 
nismaBig einfach berechnet werden kann, laBt sich diese Art der Interpolation | 
mit ertraglichem Aufwand auch fiir eine beliebige Argumentstelle durchfiihren. , 

Die Anwendung der Formel (5b**) erfordert, wenn die d,,, berechnet vor- 
liegen und [M,(w)]*+! vertafelt ist, genau so viele komplexe Multiplikationen : 
wie (5a) bzw. (5b). 

Eine Umschreibung der Formeln mit Mittelpunkt kann entsprechend vor-. 
genommen werden. 
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‘On 


Summary 


Increased attention has been recently devoted to the interpolation of an 
analytic function g(z) in the complex plane. Only if the increment in a given 
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table is sufficiently small, g(z) may be approximated by a quadratic or cubic 
polynomial. To get an interpolation polynomial P(z) of higher degree we employ 
the extension of the Lagrange formula given by HERMITE, based on the values 
of the function and some of its derivatives: 


Ms 


Lu, b (z) gl?) (Zn) < 


I 


n 
P(z) = J 
H=1 p=0 
Choosing the points z, in the corners of a regular polygon the relation for the 
coefficients L,,» may be found as follows: 


Ly,p(2) = e#(22/") P(H-1) Ty 4 fz], 
where 
glk) = g et(22/n) (1—-y), 


It is convenient and advantageous to choose a square grid. Employing e.g. 
an interpolation polynomial of the degree 11, it is unnecessary to tabulate the 
twelve coefficients L,,(z). You need only three, namely Li, Li, and Li 2. The 
corner values together with the first and second derivatives of the special function 
g(z) must be known. It should be noted, that in many cases the derivatives could 
be easily computed (i e. by the differential equation). 


(Eingegangen: 20. 4. 1951.) 
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Uber die Instabilitat von Methoden zur Integration gewéhnlicher 
Differentialgleichungen 


Von HEINz RUTISHAUSER, Ziirich?) 


Trotz der bemerkenswerten Publikation von J. Topp?), deren wesentliche 
Punkte weiter unten wiedergegeben werden, hat der Verfasser seither 6fters Ver- 
fahren zur numerischen Integration von Differentialgleichungen beobachtet, die, 
obschon zwar mit bestechend kleinem Abbrechfehler behaftet, doch die groBe 
Gefahr der numerischen Instabilitat in sich bergen. Ich will vorwegnehmen, daB 
bei den bewdhrten Verfahren von RuncE-Kutta und Apams (Extrapolations- 
verfahren), sofern sie richtig angewendet werden, diese Gefahr kaum besteht. 

Es liegt in der Natur der Sache, da8 man eine Differentialgleichung, wenn 
man sie numerisch lésen will, durch eine Differenzengleichung approximiert und 
diese lést. Um dabei die Schrittlange nicht allzuklein wahlen zu missen, bevor- 
zugt man solche Differenzengleichungen, die die Differentialgleichung méglichst 
gut approximieren, dafiir dann aber von héherer Ordnung sind als die urspriing- 
liche Differentialgleichung. Gerade darin liegt aber eine Gefahr, denn damit hat 


die Differenzengleichung auch eine gréBere Lésungsmannigfaltigkeit, und es ist 


gut méglich, daB8 die numerische Integration gerade eine der «eingeschleppten » 


1) Institut fiir angewandte Mathematik der Bs 
2) J. Topp, Solution of Differential Equations by Recurrence Relations, MTAC 4, 39-44 (1950). 
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Lésungen liefert, die nur am Anfang etwas mit der gesuchten Lésung der Diffe- 
rentialgleichung zu tun hat. In der erwahnten Arbeit von J. Topp sind mehrere 
solche Beispiele angefihrt. : 

Besonders aufschluBreich ist die Betrachtung der zugehdrigen Variations- 
gleichung: Es ist sehr gut mdglich, da6 die Differentialvariationsgleichung stabil 
ist, das hei®t nur abklingende Lésungen enthalt, wahrenddem die Differenzen- 
variationsgleichung instabil ist, indem sie dank der erweiterten L6sungsmannig- 
faltigkeit auBer den abklingenden Lésungen auch noch solche besitzt, die expo- 
nentiell ansteigen. Eine einmal vorhandene Abweichung von der richtigen Lésung, 
wie klein sie auch sein mag — und solche Abweichungen sind wegen der Aufrun- 
dungsfehler unvermeidlich —, wachst deshalb stark an und kann schlieBlich die 
gefundene Lésung betrachtlich verfalschen. Dabei ist diese Instabilitat — es sei 
hier nochmals betont — nur durch ein unzweckmaBiges Integrationsverfahren ver- 
ursacht. 

Im folgenden werden einige gebrauchliche Verfahren nach diesem Gesichtspunkt 
untersucht und einfache Kriterien fiir die Stabilitat solcher Verfahren angegeben. 
Im iibrigen aber befaBt sich diese Arbeit nicht mit Fehlerabschatzungen. 


Differentialgleichungen erster Ordnung 


i y= T(%, 9). 
Variationsgleichung: 
pen 
n ay U 


a) Integration mittels der Simpsonschen Regel?) 


h 
Vert = Vaart > (Ver + 4 Ve + Vea)*)- (1) 


Diese Beziehung liefert zusammen mit der Differentialgleichung zwei Glei- 
chungen fiir die unbekannten GréBen y;,, und y;,,,, die man meist durch Itera- 
tion lést. Ist die Differentialgleichung linear oder quadratisch in y, so kann man 
die Iteration umgehen. Wir nehmen aber an, da8 man in jedem Fall erst dann 
zum nadchsten Integrationsschritt iibergeht, wenn die Beziehung (1) erfiillt ist. 

Die Differenzenvariationsgleichung zu (1) ist offenbar 


f h 4h h 
mesa 1 ery hes) es lyk Ne — (1 ss 3 fy2-1) a= 0. 


Nimmt man /, als konstant an und macht fiir die Lésung dieser Gleichung den 
Ansatz », = 4*, so erhalt man fiir 4 eine quadratische Gleichung mit den Lé- 
sungen: ; 

h h? 


VP ong 
Wy = l+hi+ hit vos wehty, Ag = —1 + 3 fy gg fet teem —e-hfy/3 , 


: 


Dabei Sepa man leicht, da8 von den beiden Fundamentallésungen m1, = A¥ 
und 12,4 = Ay der Differenzenvariationsgleichung die erste die Lésung der Dif- 

1) Die Erscheinung wurde an diesem Beispiel auch von Herrn G. Dautoutst, Stockholm, 
beobachtet und richtig gedeutet (Vortrag an der GaMM.-Tagung 1951 in Freiburg i. Br.). Vergleiche 
auch: Z, angew. Math. Mech. 31, 239 (1951). 


i i Dabei bedeutet h die Lange des Integrationsschrittes, und es steht y;, als Abkirzung fiir 
y , 
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ferentialvariationsgleichung approximiert, wahrenddem die zweite durch das 
numerische Verfahren eingeschleppt wurde. 

Insbesondere stellt Af ~ (— 1)* e-*hiyl3 fiir f,< 0, also gerade wenn die Diffe- 
rentialgleichung stabil ist, eine langsam exponentiell anwachsende Oszillation dar. . 
Dies bewirkt, daB eine kleine St6rung der numerischen Lésung — durch einen Auf- 
rundungsfehler oder einen Abbrechfehler verursacht —im weiteren Verlauf der In- 
tegration verstarkt wird und schlieBlich vollstandig tiberhand nimmt. Die Erschei- 
nung ist bei CoLLtatz!) als Aufrauhungserscheinung bezeichnet; es werden dort 
Mittel zur Behebung des Ubelstandes angegeben. Dagegen ist die dort gegebene 
Erklarung nicht vollstandig; das Phanomen tritt nur fiir /, < 0 auf, fiir f, > 0 
ist nichts zu beftirchten, was insbesondere im Hinblick auf die gewohnliche Simp- 
sonsche Integrationsregel (f, = 0) sehr wichtig ist. 


b) Integration nach Runge-Kutta und ahnlichen Verfahren 


Da diese Verfahren y;,, aus y, nach einer eindeutigen Vorschrift und ohne 
Verwendung der vorhergehenden Werte y,_1, V,_2, -.. berechnen, bleibt die Ord- 
nung beim Ubergang von der Differentialgleichung zur Differenzengleichung un- 
verandert, es werden also keine fremden Lésungen eingeschleppt, und es ist keine 
Instabilitat zu befiirchten. 

Dieselbe Eigenschaft hat eine von W. E. MILNE angegebene Methode?). 


c) Integration nach Adams 


Wir betrachten eine Vierpunkteformel: 
h 
Yura = Vat oe (9 Vesa + 19 Yip — 5 Vea + Vea) + HP .... (2) 


Dies liefert wieder zusammen mit der Differentialgleichung zwei Gleichungen fiir 
die Unbekannten y,,, und y;,,,, die man meist durch Iteration lést. 
Die zu (2) gehérige Differenzenvariationsgleichung wird 


3h 19h 5h h 
(1 ger ae hyve) Ney — (1 re fay) "e+ a fy,e-a Mk—-1 — ore fy,k—-2 Nk-2 = 0. 


Betrachtet man wieder /, als konstant, so liefert der Ansatz n;, = A” eine Gleichung 
dritter Ordnung fiir A. Eine Lésung 4, derselben liegt sehr nahe bei e#'v, demnach 
entspricht 1, = A* der Lésung der Differentialvariationsgleichung, wahrenddem 
AE und 4% eingeschleppte Lésungen sind. 

Nun reduziert sich aber die Gleichung fiir 4 auf A? — 2? = 0, wenn h gegen 
0 strebt, so daB fiir hinreichend kleines h# jedenfalls A, und 4, klein sind, némlich 
~ + y= h f,/24 . Die eingeschleppten Lésungen 72,4 = AB und 13,4 = At klingen 
also rasch ab. Das Verfahren von ADAms ist somit fiir hinreichend kleines / stabil. 


d) Varianten von c 


Um die Genauigkeit des Adamschen Verfahrens zu steigern, werden an Stelle 
von (2) auch andere Ausdriicke fiir den Korrektor verwendet. Solange es sich 


1) L. Cotvatz, Numerische Behandlung von Differentialgleichungen (Springer, Berlin, G6ttingen 


und Heidelberg 1951). Vgl. Abschnitt 4. 4. V., Seite 51. ; i : 
2) W.E. Mitne, A Note on the Numerical Integration of Differential Equations, J. Res. nat. 


Bureau Standards 43, 5387-542 (1949). 
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dabei um die entsprechenden Fiinf- oder Sechspunkteformeln handelt, ist nichts 
dagegen einzuwenden, aber es ist Vorsicht am Platz, wenn Vz41 nicht wie bei (2) 
aus y, und den Ableitungen, sondern etwa aus y,_, oder Vz_3 und den Ableitun- 
gen berechnet wird, wie zum Beispiel bei 


2h ; ' 
Vert = Veg + 45 (7 Yieut + 32 Vig + 12 Vea + 32 Yeo + 7 Ye_s) +h" «.- 


Tatsichlich hat die zugehérige Differenzenvariationsgleichung eine Lésung 


19h 
ny = Ar mit a= (1-4 ht), 


so da8 das Verfahren fiir /,< 0 instabil ist. 


Differentialgleichungen zweiter Ordnung 


se Si taSics gh 


. Soweit diese durch Aufspaltung in ein System von zwei Gleichungen erster 
Ordnung gelést werden, gilt das bisher Gesagte. Insbesondere muB bei numerischer 
Integration von gedampften Schwingungen vor den Verfahren a und d gewarnt 
werden. 

Es gibt aber auch Verfahren, welche eine Gleichung zweiter Ordnung ohne 
Verwandlung in ein System auflésen: 


e) Das Verfahren der zentralen Differenzen) 


Die Formeln, die diesem Verfahren zugrunde liegen, sind (speziell fiir zweite 
Ordnung): 


. h? u ” ” 
Vest = 2Ie— Vert Ty (View + 10 Ye + Ver)» (4) 
t , h v ” ” 
View = Vena + > (Yeu + 496 + Yea) - (5) 


Sie liefern zusammen mit der Differentialgleichung drei Gleichungen fiir die Un- 
bekannten Viva» Vex und yx,,. Die zu (4) und (5) gehérigen zwei simultanen 
Differenzenvariationsgleichungen werden unter der Annahme von konstantem f, 
und fy mit dem Ansatz n, = p A*, ), = q A* gelést; man erhalt wegen 

Neva = fy Nee t+ fy News 
mit den Abkiirzungen a fiir h® f,/12 und b fiir h f,/3 die Gleichungen 


p((A® — 20+) — a (42 4102+ 1) — g 4b (a8 + 10. +1) = 0 [aus (4)] 


(6) 
q [(A? — 1) —b(M 44441) — pia (A+ 4441) = 0 [aus(5)). 


Diese kénnen mit (p, q) + (0,0) nur dann simultan bestehen, wenn die Deter- 
minante dieses Gleichungssystems fiir p und g verschwindet; man erhalt nach 


1) Vergleiche L. Cotratz, Numerische Behandlung von Differentialgleichungen, S. 80. 
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einiger Rechnung: 

= 1) (A? — 1) (A— 1) — a (A+ 1) (A424 1044-1) —b(A— 1) (A?+4A4+41)])=0. 
Die vier Lésungen dieser Gleichungen sind 


A= 1lt+ah+-::: wobei a, und a, die Lésungen der Gleichung 
Ag=1+ah+--- J ot?§—afy —f, = 0sind, 
Tes ik 


“eset Seat 


Offenbar sind AF und AF die regularen Lésungen der Differenzenvariationsglei- 
chung; sie entsprechen zwei Fundamentallésungen der Differentialvariationsglei- 
chung, dagegen sind A; und A, eingeschleppt. Solange fy = 0 ist, ist nichts zu be- 
fiirchten, insbesondere kann man das Verfahren fiir eine y’-freie Gleichung sehr 
empfehlen, aber fiir fy < 0 (gedampfte Schwingungen) wdchst 


AS 
I 


yp Aes (= 1)® e-(kh/3) ty 


an, und auch A, kann noch gefahrlich werden, weil auch A3, = 1 relativ zu einer 
gegen Null konvergierenden Funktion schlieBlich sehr groB wird. 
Der Verfasser hat das Beispiel y” + y’+ 1-25 y =0 mit den Anfangsbedingungen 
Yo = 0, vp = 1 (exakte Lésung: e-+/2 sin x) auf der programmgesteuerten Rechenmaschine 
der ETH. durchgerechnet. 
Es folgen einige Ausschnitte aus der so erhaltenen Funktionstabelle (es wurde mit 
kh =0,1 gerechnet): 


, 


% y, ay 
4,8 — 0,0903699 | 0,0531227 
4,9 — 0,0847792 | 0,0584842 
5,0 —0,0787132 | 0,0626410 
Syl — 0,072 2891 0,065 6573 
D2 — 0,0656173 0,067 6070 


In dieser Region ist noch nichts Auffalliges zu bemerken, die y-Werte weichen ungefahr 
um eine Einheit der letzten Stelle von den wahren Werten ab, und nur Differenzenbildung 
bei den y’-Werten enthiillt eine gewisse UnregelmaBigkeit. Bei ¢ = 17 dagegen zeigt sich 
der Einflu8 von A bei den y’-Werten schon stark und auch bei den Differenzen der y- 


Werte: 


z y | y’ 
17,0 — 0,00019574 | 0,00005017 
ibis — 0,00019061 | 0,00005253 
L752 — 0,00018366 | 0,000086 20 
73 — 0,000175 24 | 0,00008239 
17,4 — 0,00016548 | 0,00011235 
ED — 0,00015475 | 0,00010258 | 


|, i a ee ee SS ee 
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Das wesentlich schwichere Oszillieren der y-Werte folgt auch aus den Gleichungen (6); 
fiir A = J, erhalt man namlich aus der ersten dieser Gleichungen Ww 


p ae =e wets 
aa 45> fy’, hier also p 600° 


Der weitere Verlauf der numerischen Integration bedarf keines Kommentars: 


—0,00000815 0,000 053 20 
—0,00000864 | —0,00006078 
—0,000008 62 0,000 05968 
—0,00000887 | —0,00006247 
—0,00000861 0,000 066 01 
—0,00000868 | —0,00006486 


—0,00000140 | —0,00053037 


0,000 00041 0,000 54889 ; 
—0,00000144 | —0,00056693 y 
0,000 00050 0,000 58682 
—0,00000148 | —0,00060603 
0,000 00060 0,000 627 35 


Der Verfasser ist sich wohl bewu8t, daB die Annahme von konstantem /, und 
fv in den obigen Betrachtungen die Allgemeinheit wesentlich einschrankt. In- 
dessen zeigen die Ergebnisse, daB es anscheinend nur auf das Vorzeichen dieser 
Gr6Ben ankommt, und dieses ist in der Tat in sehr vielen Fallen unverdnderlich. 
Die Aussagen sind deshalb qualitativ fast allgemeingiiltig. Nur wenn zum Bei- 
spiel bei Methode a /, im Verlauf der Integration sein Vorzeichen von Zeit zu_ 
Zeit andert, tritt ein besonderer Fall ein, indem die auftretenden Oszillationen 
abwechselnd anschwellen und wieder abklingen. 


Allgemeine Betrachtungen 


Die Betrachtung der Beispiele legt die Vermutung nahe, da8& Instabilitat 
genau bei solchen Integrationsverfahren auftreten kann, welche 4; 4, durch Inte- 
gration von y’ tiber mehrere Intervalle hinweg bilden (zwei Intervalle bei Ver-— 
fahren a, ein Intervall bei b und c, 4 Intervalle bei d, 2 Intervalle bei e). Dies 
trifft jedoch nicht genau zu, und wir wollen deshalb eine allgemeine Integrations- 
methode einer Priifung unterziehen: 

Fast alle bekannten Verfahren beniitzen Beziehungen, welche in der allge- 
meinen Formel enthalten sind: 
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(0) (0) ay" N 
Pia = 2 Gj Yeap oe hs Ws Dy gr oo eh > bay Vins 
= —m 


Ve = » al) Vaagt h Ss ay) Ving t hs te a an) VEY 


RE Bo memo me mnL LSS 0 10 ke) (ene (= 1s) 14) (6) 610. [erse ese) 161. 0).6) Welfare (6° ~ (0) ei (| 6) “ele sil) 18\ (6) 0 (6 e/e, 0) wie) is eriel 6 eleven 


Dazu kommen Differentialgleichungen 


(1) te 
eee SURO co Vavs) = 0), 
(n+1) 
Fite Verve so Views )= vs (8) 
(N) aS 
Find ered yates Vea een 


welche zusammen mit den Gleichungen (7) N + 1 Gleichungen fiir die N + 1 Un- 


bekannten Yyi1, View oo ye bilden. Im allgemeinen ist N = n, doch beniitzt 
W. E. MILNE in der erwahnten Methode!) héhere Ableitungen, als in der Diffe- 
rentialgleichung auftreten. Er differenziert diese deshalb mehrmals, um die 
notige Anzahl von Beziehungen zu erhalten. Man kann also die Gleichungen 
Fy4,+-. bis Fy durch Differenzieren der Ausgangsgleichung F,, erhalten. 

Wenn man in den Gleichungen (7) noch alles auf eine Seite bringt, so lauten die 


Variationsgleichungen fiir das ganze System (7) und (8) offenbar (mit at) =a 
yl . . 
x FY all pe- fh) =—Or (¢ = 0,1, ..., % — 1) 
as “wy 
w=ij=—m 
4 > 
Fi nw ails (j=n,n+1,..., N) 
0 oy) 


Macht man hiefiir den Ansatz nf) = Pp, #', so erhalt man beim Einsetzen ein 
System von N + 1 simultanen Gleietangen fiir die p,, die nur erfiillt sein kén- 
nen, wenn die Determinante des Systems verschwindet: 


N , 
ses Se 97) Wee? = 0, (6 sey Le, aa) 
pat \G=—m 

i ; 

Ee == (5 (a el aera N) 


Definiert man noch mit den in den Formeln (7) auftretenden Koeffizienten ar ? die 
Funktionen 


1) FuBnote 2 auf Seite 67. 
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wobei Ai, = 0 fiir uw <7, so lautet die charakteristische Determinante wie folgt: 


Aber WAg Aga Wee cies oe os ene rae hN Aon 

0 EE Jewelry cocoon Sneak. 5s aucomonons . AN-1 Ain 

0 0 0 An-1,n-1 hN-"+1 An-1,N 
DA) we | ae oie Se ata les is btn aren siaimie mw © om ee <n en slo Sina g a oe ee 

Fn, y Fn, 9’ Fy, y™ 0 0 

Fn+ij,y Fr+1,9’ : 

Fy,y Fy, Y", areiie Se Spe eis o)iaue, oo Rim es # us legate’ silo) @ ie ieale arene ji Fy, y(N) 


Wenn das Verfahren ein Polynom i-ten Grades, das die Lésung einer Differen- 
tialgleichung ist, samt seinen Ableitungen exakt wiedergeben soll, und dies darf 
man fordern, so miissen die Bedingungen 


. : pM oy eee a ye 
yi) = : yd) ye Lynas abl = yf J=0 


vertraglich sein. Daraus folgt aber durch Einsetzen in die i-te der Gleichungen 
1 


Die Gleichung D(A) = 0, die iiber die Stabilitat des Verfahrens entscheidet, 
mu8 nun ~ Lésungen in der Nahe von A = 1 haben, entsprechend den  unabhan- 
gigen Lésungen der Variationsgleichung. In der Tat findet man fiir h = 0, wo 
sich D(A) bis auf einen Faktor auf Ao9* 411 +++ An-1,n-1Teduziert, daB A = 1 wegen 
A,,; (1) = 0 eine n-fache Nullstelle von D(A) ist. 

Alle andern Nullstellen von D(A) entsprechen damit eingeschleppten Lésungen 
der Differenzengleichung und miissen, damit das Verfahren stabil ist, fiir hin- 
reichend kleines h im Innern oder héchstens auf dem Rand des Einheitskreises 
liegen. Dies ist sicher dann der Fall, wenn fiir h = 0 alle Nullstellen von D(A) im 
Innern des Einheitskreises liegen, und sicher nicht, wenn einzelne auBerhalb | 
sind. Somit 

Hinveichende Bedingung fiir Stabilitaét des Verfahrens (7) fiir geniigend kleines h: 

Samtliche Funktionen 

A (a) = sy ait) qitm 


besitzen auBer der trivialen einfachen Nullstelle 4 = 1 nur noch solche Nullstellen 
mit As) als 

Notwendige Bedingung: Keine der Funktionen A,;,(A) besitzt eine Nullstelle 
auBerhalb des Einheitskreises?). 


Die tubrigen Funktionen A;,(A) kénnen die Stabilitat nur beeinflussen, wenn 
zwar die notwendige, aber nicht die hinreichende Bedingung erfiillt ist. 


1) J. Topp hat Verfahren betrachtet, welche nicht einmal die notwendige Bedingung erfillen. 
Die erhaltene Loésung wird dann schon nach wenigen Schritten véllig unbrauchbar. 
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Anwendungen 


Fur die Formel (1) ergibt sich (es ist N= n= 1,m= 1): 


il 
Ang = At 4-1, Borie AS EA Al)”. = 9-H, 9) > 


Somit 
h 
det A2y 7 (A2 ad 1 
D(a) = 3 Aaa Aaa) 
—f, 1 
DaB Ao zwei Nullstellen mit |4| = 1 hat, mahnt schon zur Vorsicht, aber auBer- 
dem liest man sofort ab, daB D(A) fiir {,< 0 und A = —1 positiv ist, fiir A = — oo 


dagegen negativ. Somit liegt eine Nullstelle links von —1; das Verfahren ist 


instabil. 
Fiir die Formel 5.42 im erwahnten Buche von CoLramz (S. 81) ist (V = = 4, 


m = 1) 
FINA IS 9%, 


Alle andern A;,, kommen nur mit mindestens h? in der Determinante vor. Somit 
ist bis auf Glieder mit h?: 


(A ae 0 0 0 0 
0 —APF + 1 2Ah 0 0 
D(a = 0 0 eae 0 0 
0 0 0 =A? 1 Fats 4a41) 
—fy =/,' ae? —f,/” i 
Fiir 2 = — oo ist D positiv, fiir 1 = —1— ¢« hat D das Vorzeichen von 


—A2set = (A2 + 42 +1) 
—fy” it 


welches fiir fy < 0 und hinreichend kleines ¢ negativ ist. Also ist dieses Verfahren 


-instabil fiir fy”< 0. mei 
Anderseits ist es leicht, Verfahren anzugeben, die immer stabil sind. Man 


braucht nur die Formeln (7) so zu gestalten, daB jede Zeile mit 


t+1 


1 
ro(4 ; i ( i =5 
Dates Vat iS Alves Say las (i = 0, 1,...,%— 1) 


beginnt. Damit ist A,,(4) = — A"+1+ 4™ und hat somit nur die triviale Nullstelle 
4 = 1 auf dem Rande des Einheitskreises. 
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Summary 


In the numerical solution of a differential equation as a difference equation, 
the latter is usually of higher order and has therefore more solutions than the ori- 
ginal differential equation. It may well be that some of these “‘extra”’ solutions 
grow faster than any solution of the given equation; in this case the computa- 
tional solution has the tendency to follow one of these and has after a certain 
number of integration steps nothing to do with the original differential equation. 

The author gives some examples and a criterion for stability of integration 
methods. This criterion is then applied to some well-known integration formulas. 


(Eingegangen: 25. 9. 1951.) 
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8. Internationaler Kongre® fiir theoretische und angewandte Mechanik 


Der Kongre8 wird vom 20. bis 28. August 1952 an der Universitat Istanbul 
stattfinden und die Sektionen 
I. Elasticity, Plasticity, Rheology, 
II. Fluid Mechanics (Aerodynamics, Hydrodynamics), 
III. Mechanics of Solids (Ballistics, Vibrations, Friction, Lubrication), 
IV. Statistical Mechanics, Thermodynamics, Heat-Transfer, 
V. Mathematics of Physics and Mechanics, Methods of Computation 
umfassen. Interessenten erhalten Auskunft (insbesondere auch iiber die vorge- 
sehenen Reiseerleichterungen) durch: The Eighth International Congress for 
Theoretical and Applied Mechanics, P. O. Box 245, Istanbul, Turkey. 
H. Favre und H. ZIEGLER 


4. KongreB der IVBH. 1952 


Die Internationale Vereinigung fiir Briickenbau und Hochbau fihrt ihren 
4. KongreB vom 24. August bis 1. September 1952 in Cambridge und London 


durch. Es werden sechs Arbeitssitzungen durchgefiihrt, die den folgenden Themen 
gewidmet sind: 


A. Allgemeine Fragen. 

I. Bemessungsgrundlagen und Sicherheit. 

II. Entwicklung der Berechnungsmethoden. 
B. Stahlbau. 

I. Grundlagen. 

II. Praktische Anwendungen, 
C. Massivbau. 

I. Grundlagen und Eigenschaften des Betons. 

II. Aktuelle Probleme des Betons und des Eisenbetons; vorgespannter Beton. 

Teilnahme am KongreB setzt die Mitgliedschaft der IVBH. voraus ( Jahres- 
beitrag Fr. 10.— fiir Einzelmitglieder) ; Anmeldungen sind zu richten an das 
Sekretariat der Internationalen Vereinigung fiir Briickenbau und Hochbau 
(IVBH. 1952), ETH., Ziirich. F. Stuss1 


—— 
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Response of Physical Systems. By JouHNn D. Trimmer (John Wiley & Sons, 
New York 1950). 268 pp., 93 figs.; $5.—. 


_ Das Buch beginnt mit einer sehr allgemein gehaltenen Diskussion der folgen- 
den Art. Ein beliebiges System sei irgendeiner Einwirkung (forcing) ausgesetzt 
und erzeuge eine bestimmte Wirkung (response). Als Beispiel zitiert der Ver- 
fasser die Familie Jones, die von einer Stadt nach einer anderen umzieht (forcing) 
und durch diese Einwirkung ihre Sprachgewohnheiten 4ndert (response). Ein 
_solch allgemeines Schema, in das selbst die Anderungen der Sprachgewohnheiten 
der’ Familie Jones hineinpassen, 1a8t sich ohne Zweifel aufstellen, doch erscheint 
es fraglich, ob dadurch unsere Erkenntnis in irgendeiner Form erweitert wird. 
Auch der Verfasser mu8 nach dieser sehr allgemeinen Einfiihrung zu «iiblichen», 
wohldefinierten Systemen iibergehen. In der Folge werden eine Reihe verhdltnis- 
maBig einfacher Anordnungen aus verschiedenen Gebieten der Physik auf ihr 
Verhalten untersucht, falls wir eine Wirkung auf ein solches System ausiiben. 
Die sich’gleichbleibenden Methoden der mathematischen Behandlung bilden die 
organische Verbindung zwischen den einzelnen diskutierten Gebieten. 
Mit einem kurzen Abschnitt tiber nichtlineare Systeme und tiber die Verwen- 
dung der Laplace-Transformation zur Behandlung linearer Anordnungen schlieBt 
das Buch. E, Baldinger 


Vorlesungen tiber héhere Mathematik. Bd. 2. Von A. DUSCHEK (Springer- 
NWerlag, Wien 1950). 386 S.; sFr. 33.90. 

In diesem zweiten Band der Vorlesungen DuscHEKs sind die unendlichen 
Reihen, Funktionen mehrerer Veranderlichen, einige fundamentale Ergebnisse 
der Wabhrscheinlichkeitsrechnung und ihre Anwendung auf Fehlertheorie und 
Ausgleichsrechnung sowie die lineare Algebra und Tensorfelder dargestellt. Die 
Vorlesungen sind im wesentlichen fiir solche Zuhérer bestimmt, die Mathematik 
nicht als Selbstzweck betreiben, wie Ingenieure, Physiker, Chemiker und andere. 
Aus diesem Grunde werden die Resultate vor der Darstellung der exakten Beweise 
zuerst nach Moéglichkeit anschaulich begriindet. Aufgaben samt Lésungen er- 
lauben es dem Leser, sich am dargestellten Stoff zu tiben. 

Die Darstellung des Stoffes ist sauber und klar, seine Auswahl fiir solche Vor- 
lesungen von allgemein orientierendem Charakter zweckmaBig. Auch Ergebnisse 
der modernen mathematischen Forschung wurden geschickt verarbeitet. Aus 


allen diesen Griinden diirfen diese Vorlesungen bestens empfohlen werden. 
W. Saxer 


Praktische Funktionenlehre, Bd. 1: Elementavre und elementare tvanszen- 
dente Funktionen. Von F. TOLkeE (Springer-Verlag, Berlin 1950). 440 S., 178 Abb.; 
DM 39.-. 

Der erste Teil des Bandes befaBt sich mit den sogenannten elementaren Funk- 
‘tionen; es werden Differentialgleichungen, denen diese Funktionen gentigen, und 
andere vielgebrauchte Beziehungen zusammengestellt. Ferner enthalt dieser Ab- 
schnitt eine reiche Auswahl an unbestimmten Integralen von trigonometrischen, 
hyperbolischen, zyklometrischen und Wurzelausdriicken, welche dem Praktiker 
sehr willkommen sein diirften. Den Abschlu8 des ersten Teils bilden Tabellen 


der elementaren Funktionen. 
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Im zweiten Teil wird gezeigt, wie man unter Herbeiziehung von Fourier- 
Reihen ausgewahlte Aufgaben der Technik lést. Soweit bei diesen Aufgaben neue 
Funktionen auftreten, werden sie eingehend behandelt und tabelliert. Insbeson- 
dere werden Thetafunktionen und ihre Integrale durch Betrachtung von Warme- 
leitungsaufgaben und die Kugelfunktionen an Hand von Aufgaben tiber bela 
Membrane eingefiihrt. Ferner werden Platten- und Balkenprobleme behandelt. 

Neuartig an diesem Werk ist vor allem die durch Vereinigung von Theorie, 
Beispielen und Funktionstafeln erreichte Anschaulichkeit. H. Rutishauser 


Vorlesungen iiber Differential- und Integralrechnung. Von L. VIETORIS 
und G. Locus (Universitatsverlag Wagner, Innsbruck 1951). 415 S., 159 Abb.; 
sFr. 28.-. 

Das Buch vermittelt, unter Beschrankung auf das Wesentliche, eine Einfiih- 
rung in die Differentialrechnung der Funktionen mit einer und mehreren Varia- 
blen und in die Integralrechnung der Funktionen mit einer Variablen, umfaBt also 
den Stoff, der normalerweise in den beiden ersten Semestern geboten wird. Einge- 
schaltet ist noch ein Kapitel iiber Vektorrechnung. Einige Ubungsaufgaben 
dienen der Festigung des Wissens. Bemerkenswert ist die Sorgfalt, mit der neue 
Begriffe eingefiihrt werden und die Beispiele durchgerechnet sind. Eine gewisse 
Breite der Darstellung erleichtert dem Studierenden — an den sich das Buch, als 
sehr empfehlenswerte Hilfe, wendet — das Erarbeiten des Stoffes. 

E. Roth-Desmeules 


Der Kreisel; seine Theorie und seine Anwendungen. Von R. GRAm- 
MEL (Springer-Verlag, Berlin 1950). 1. Band: Die Theorie des Kreisels, 281 S., 
137 Abb.; 2. Band: Die Anwendungen des Kreisels, 268 S., 133 Abb.; DM 30.— 
pro Band. 


In der umfangreichen Literatur auf dem Gebiete der Mechanik ist der Kreisel 
(und mit ihm die héhere Kinetik des starren Kérpers iiberhaupt) von jeher mit 
nur wenigen Werken vertreten gewesen. Um die Jahrhundertwende schufen 
F. KLEIN und A. SOMMERFELD mit ihrer Theorie des Kreisels das heute noch 
anerkannte, aber sehr abstrakt gehaltene Hauptwerk. Im Jahre 1920 stellte ihm 
R. GRAMMEL mit seinem Kveisel ein Buch zur Seite, das sich mit seiner betont 
anschaulichen Haltung und seiner Darstellung der wichtigsten technischen An- 
wendungen an den Ingenieur richtete. Dieses Werk, das groBer Nachfrage begeg- 
nete und schon seit Mitte der dreiBiger Jahre vergriffen ist, liegt nun in einer 
zweibandigen, stark umgearbeiteten und modernisierten Neuauflage vor. : 

In der neuen Gliederung umfa&t der erste Band im wesentlichen die Theorie, — 
der zweite die technischen Anwendungen des Kreisels. Im theoretischen Band 
wird nach einer ausfiihrlichen Grundlegung aus didaktischen Griinden erst der 
symmetrische und dann der unsymmetrische Kreisel behandelt. Ein weiteres 
Kapitel befaBt sich mit allgemeineren Problemen und enthalt insbesondere eine 
sehr wertvolle grundsatzliche Untersuchung gyroskopischer Systeme. Ein Anhang 
ist der Darstellung der Bewegung durch elliptische Funktionen gewidmet. Der 
zweite Band gibt eine reichhaltige Ubersicht der technischen Probleme und An- 
wendungen. Er behandelt die (meist unbeabsichtigten) Kreiselwirkungen bei — 
Radsatzen, wie sie insbesondere im Fahrzeug- und Flugzeugbau vorkommen, 
sodann die Kreiselgerate (Kreiselkompasse, kiinstliche Horizonte, Wende- und 
Lagekreisel) und schlieBlich die unmittelbaren gyroskopischen Stabilisatoren. 

Klar, anschaulich und doch kompromiflos streng geschrieben, im Theoreti- 
schen bis zu den bekannten Lésungen fiihrend und ebenso griindlich in der 
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Behandlung der technischen Probleme, diirfte dieses Werk auf langere Zeit das 
Lehrbuch vom Kreisel werden. Ai, Ziegler 


Methoden der praktischen Analysis. Von F. A. WILLERS (Walter de 


_ Gruyter & Co., Berlin 1950) 410 S., DM. 24.-. 


Das vorliegende Buch stellt eine verbesserte und erweiterte Auflage des 
1928 zum erstenmal erschienenen Werkes dar. Neu sind die Abschnitte iiber die 
Integration von Integralgleichungen zweiter Art, tiber die Auswertung uneigent- 
licher Integrale sowie tiber die iterative Lésung linearer Gleichungssysteme. 
Dariber hinaus wurde das ganze Kapitel VI mehr den heutigen Bediirfnissen an- 
gepaBt. Im ersten Kapitel wird auf die Grundlagen und Hilfsmittel des Zahlen- 
rechnens eingegangen, wobei auch die Nomographie zu Worte kommt. Mit beson- 
derer Sorgfalt werden im zweiten Kapitel die Interpolation und die damit zusam- 
menhangende numerische Integration und Differentiation behandelt. Sodann 
folgen Betrachtungen iiber graphische Integration unter Einschlu8 mechanischer 
Hilfsmittel. Kapitel IV ist der Ausgleichsrechnung (Methode der kleinsten Qua- 
drate) und der harmonischen Analyse gewidmet, Kapitel V den héheren Glei- 


-chungen und linearen Gleichungssystemen. In Kapitel VI schlieBlich sind die 


Verfahren zur genaherten Integration von gew6hnlichen Differentialgleichungen — 
insbesondere die Methoden von RuncE-Kutra und Apams — dargestellt, wobei 
auch partielle Differentialgleichungen und Eigenwertprobleme gestreift werden. 

A, Rutishauser 


Anwendungen der elliptischen Funktionen in Physik und Technik. 
Von F. OBERHETTINGER und W. MaaNnus (Springer-Verlag, Berlin 1949). 126 S., 


@54 Abb., DM. 15.60. 


Der Akzent liegt auf Anwendung. Auch auf eine kurze Einfiihrung in die 
Theorie konnten die Verfasser dank der vorhandenen ausgezeichneten Lehrbticher 
verzichten. Das erste Kapitel bringt deshalb lediglich eine knappe Zusammen- 
stellung der Definitionen und Formeln elliptischer Funktionen nach JAKoBI und 
WEIERSTRASS, der Thetafunktionen und der elliptischen Integrale. Das zweite 
Kapitel behandelt die konformen Abbildungen, die den elliptischen Funktionen 
zuganglich sind, vor allem polygonaler Gebiete, und die Berechnung der zugehdri- 
gen Greenschen Funktion. Die Anwendungen auf elektrostatische Probleme 


(3. Kapitel) betreffen vor allem die Kapazitaét zweidimensionaler Kondensatoren 
und das Feld von linearen und quadratischen Drahtgittern. Die aerodynamischen 


Beispiele (4. Kapitel) handeln von der Bewegung von Wirbelsystemen und vom 
Tragfliigel im Windkanal, wahrend das letzte Kapitel Beispiele vermischten In- 
haltes bringt, von denen die Tschebyscheffsche Approximation in einem fiir die 
Theorie der elektrischen Filter wichtigen Fall besonders zu erwahnen ist. 

Die verschiedenen Fragen sind sehr sorgfaltig behandelt und bis zum formel- 
maBigen SchluBresultat durchgefiihrt. So enthalt das Buch eine Menge von An- 
regungen wie durchgerechneter Beispiele. A, Pfluger 


Vorstufe zur theoretischen Physik. Von R. BEecKER (Springer-Verlag, 
Berlin 1950), 172 S., 94 Abb.; DM 7.50. 

Die erste Bekanntschaft mit einem neuen Buch vollzieht sich meist in der 
Weise, da8 man sich das Inhaltsverzeichnis ansieht. Im vorliegenden Fall erlebt 
man dabei gleich eine Uberraschung. Die vier Hauptabschnitte des Buches tragen 
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: 
die Uberschriften ; Aus der Mechanik; Schwingungen und Wellen; Aus der Warme- 
lehre; Mathematische Erinnerungen und Beispiele. Wo, wird sich der Leser 
fragen, bleiben die Elektrizitatslehre, die Optik, die Atomphysik usw. ? Was will 
das Buch, in dem einige willkiirlich herausgegriffene spezielle Probleme der Me- 
chanik und der WArmelehre ausfiihrlich besprochen und durchgerechnet sind ? 
Offensichtlich will es kein vollstandiges Lehrbuch der theoretischen Physik sein, 
von denen es viele gibt; der Verfasser verfolgt vielmehr ein anderes Ziel. Es geht 
ihm darum, bei jungen Physikstudenten Sinn und Verstandnis fiir die mathe- 
matische Formulierung physikalischer Tatbestande und Gesetze zu wecken. Die 
Notwendigkeit eines Buches, welches diesem Zwecke dient, ergibt sich aus der. 
Tatsache, daB einerseits in mathematischen Vorlesungen und Ubungen oft zu 
wenig Gelegenheit geboten wird, konkrete physikalische Probleme zu diskutieren, 
wahrend andererseits in den Vorlesungen iiber Experimentalphysik ihre strenge 
mathematische Behandlung zu kurz kommt. Uber diese je nach Hochschule mehr 
oder weniger groBe Kluft schlagt das Buch eine Briicke; und dies geschieht mit 
Erfolg. . 
Die behandelten Beispiele sind gut gewahlt. Der Text ist klar und auBerst 
lebendig abgefaBt, und der Leser wird in anregender Art pers6nlich angesprochen 
und zu Mitarbeit und Nachdenken aufgefordert. Jeder Student, der den darge- 
botenen Stoff griindlich durcharbeitet, wird groBen Nutzen daraus ziehen und 
sich den Zugang zur reinen theoretischen Physik erheblich erleichtern. G. Busch 


7 
; 


Die Welt der Vektoren. Von F. OLLENDORFF (Springer-Verlag, Wien 1950), 
470 S., 68 Abb.; sFr. 39.-. 


Wie der Verfasser im Untertitel anfiihrt, umfaBt das vorliegende Werk eine 
Einfitihrung in die Theorie und Anwendung der Vektoren, Tensoren und Opera- 
toren. Es ist die Frucht der Erkenntnis, daB fiir das Studium der exakten Natur-_ 
wissenschaften, insbesondere aber der theoretischen Physik, die bisher iibliche 
Basis der Vektorlehre zu schmal geworden ist. InhaltsmaBig nimmt das Werk 
eine Mittelstellung ein zwischen einem Lehrbuch der reinen Mathematik und 
einem Lehrbuch der theoretischen Physik. 

Die einzelnen Kapitel befassen sich mit folgenden Gegenstanden: Darstellung 
der elementaren Vektorrechnung im dreidimensionalen Raum; Vektoranalysis. 
tiber der euklidischen Gruppe; Vektoralgebra im affinen Raum von » Dimen- 
sionen; Tensoralgebra und Tensoranalysis iiber der affinen Gruppe; Darstellung 
der Grundlagen der Riemannschen Geometrie; Einfiihrung in die spezielle Rela- 
tivitatstheorie; der Hilbertsche Raum und die Theorie der Operatoren (hier 
quadratische komplexgliedrige Matrizen). In die einzelnen Kapitel sind zahl- | 
reiche langere und kiirzere Abschnitte iiber Anwendungen des Vektor- und 
Tensorkalkiils in der theoretischen Physik, speziell der Elektrodynamik und der 
Elastizitatstheorie, eingeflochten. 

Das Buch ist fiir fortgeschrittene Leser geschrieben; bereits im ersten Kapitel 
wird ohne weitere Erklarungen die Sprache der Indizes verwendet. Durch die 
Zusammenfassung von Theorie und Anwendung des Vektor- und Tensorkalkiils 
verschafft das Werk einen guten Uberblick iiber diese mathematische Disziplin, 
welche zum Verstandnis des heutigen physikalischen Weltbildes unbedingt not- 
wendig ist. M. Jeger 
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herausgegeben v. Dr. HERMANN EBERT 
(Physikalisch-Technische Bundesanstalt) 
1951. KI. 8°. 528 S. Diinndruck mit 147 Abb., 
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Ein handliches Nachschlagewerk, wie es 
bisher gefehlt hat. Mit seinen kurzen, pra- 
gnanten und doch umfassenden Begriffs- 
bestimmungen, Zahlenwerten und MeBda- 


ten gibt es alle Grundlagen fiir wissen- 
schaftliches Arbeiten. Dem in Wissenschaft 
und Praxis arbeitenden Physiker ebenso wie 
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Taschenbuch eine ganze Bibliothek. 
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Mitteilungen aus dem Institut fiir angewandte Mathematik der ETH., Nr.1 und 2 


Dr. A. P. SpeIser, ETH., Ziirich 
Entwurf eines elektronischen Rechengerates 


Unter besonderer Berticksichtigung der Erfordernis eines minimalen Material- 
aufwandes bei gegebener mathematischer Leistungsfahigkett. - 70 Seiten mit 3 Fi- 
guren. Kartoniert Fr. 6.75 (1950). 

Das Buch ist aus dem Bediirfnis heraus entstanden, einen programmgesteuerten 
Rechenautomaten fiir die Verwendung in einem kleinen bis mittleren mathema- 
tischen Institut zu schaffen. 


Dr. H. RutisHAusEr, Dr. A. P. SpeIser, Prof. Dr. E. Stieret, ETH., Ziirich 


Programmegesteuerte digitale Rechengerate 


102 Seiten mit 28 Figuren. Kartoniert Fr. 8.85 (1951). — Inhalt: Grundlagen und 
wissenschaftliche Bedeutung — Organisation und Arbeitsweise — Arithmetische Prin- 
zipien — Vorbereitung von Rechenplanen — Physikalische Grundlagen. 

Damit wird endlich eine empfindliche Liicke im deutschsprachigen Schrifttum 
geschlossen, das bisher auBer Referaten und Zeitungsnotizen keine wissen- 
schaftliche Zusammenfassung der Entwicklung von GroBrechenautomaten ~ 
aufzuweisen hatte. 
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Neuerscheinungen 


Theoretische Elektrotechnik 


Band III: Grundztige der Theorie elektrischer Maschinen. Von Dr. Ing. Karu 
KUHLMANN, em. Professor an der Eidgenéssischen Technischen Hochschule in 
Zirich. (Lehr- und Handbiicher der Ingenieurwissenschaften, Band 13.) 546 Sei- 
ten mit 328 Abbildungen. In Ganzleinen Fr. 74.90. 


Elektrische Maschinen 


Band I: Allgemeine Berechnungselemente. Die Gleichstrommaschinen. Von Dr. 
Rupo.r Ricuter, em. Professor an der Technischen Hochschule in Karlsruhe, 
vorm. Direktor des Elektrotechnischen Instituts. 

Zweite, verbesserte Auflage 1951. 642 Seiten mit 4538 Textabbildungen. In Ganz- 


leinen Fr. 48.-. 
Elemente der Philosophie und der Mathematik 


Eine Anleitung zum inhaltlichen Denken. Von Dr. ANDREAS SPEISER, Professor 


an der Universitat Basel. (Wissenschaft und Kultur, Band 6.) 116 Seiten. In 
Ganzleinen Fr. 11.45, 
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